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Resumen 


El grafeno es un alótropo bidimensional del carbono, con propiedades físicas y 
químicas excepcionales, es un material flexible, liviano, transparente, con alta 
movilidad de carga y gran área superficial, entre otras. Desde su obtención ex- 
perimental en el 2004 por Andre Geim y Konstantin Novoselov, hasta la fecha, 
ha sido un material prometedor con múltiples aplicaciones en campos como la 


electrónica, bioingeniería, biomedicina y nanotecnología. 


Particularmente, en la construcción de biosensores, el grafeno es un material que 
ofrece muchas ventajas, ya que, es sensible, selectivo, de bajo ruido electróni- 
co intrínseco y además, biocompatible. Trabajos experimentales y teóricos dan 
cuenta de la facilidad del grafeno para adsorber bacterias, células e incluso 
moléculas, hecho que permite estudiar qué propiedades específicas modifica 
el material luego de la adsorción, información útil para la eventual funciona- 
lización y caracterización de un biosensor basado en grafeno. Trabajos en esta 
dirección, con nucleobases de ADN han sido explorados por diversos investiga- 


dores. 


Dependiendo del tipo de interacción generada entre el grafeno y las moléculas 
de interés biológico, los mecanismos o modos de sensado pueden ser: ópticos, 
magnéticos, mecánicos y electrónicos. De acuerdo a la literatura, en los últimos 
años ha ganado gran interés el desarrollo de transistores de efecto de campo ba- 
sados en grafeno (GFETSs), que son dispositivos electrónicos con tres contactos 
metálicos (fuente “source”, drenaje “drain d” y compuerta “gate”) y un canal de 
grafeno. Al aplicar voltajes de polarización tanto en el drenaje como en la com- 
puerta, es posible medir la respuesta del dispositivo, obteniendo la conductancia 
y las curvas de corriente características, antes y después de la interacción entre 


la superficie de grafeno con biomoléculas. 


En este trabajo se estudiaron teóricamente los efectos que produce la fisisorción 
de diez aminoácidos en las propiedades electrónicas del grafeno, para evaluar 
su eventual aplicación en un detector de aminoácidos. Para ello se modeló un 
dispositivo que comprende dos electrodos de grafeno (fuente R y drenaje L) 
conectados a un canal de grafeno. Dependiendo del aminoácido adsorbido sobre 
el canal, se pueden modelar los efectos en las propiedades eléctricas, debido a 
las diferencias en la transferencia de carga, la concentración de dopaje, el tipo de 
dopante y la generación de dipolos eléctricos locales. Se calcularon las curvas de 
transmisión y de corriente-voltaje características de cada sistema, considerando 


voltajes de polarización positivos y negativos (entre -2 y 2 V). 


Esta tesis se dividió en seis capítulos a saber: 1) Introducción, 11) Conceptos 
teóricos, 111) Descripción de los sistemas modelados, 1v) Resultados y discusión 
de adsorción de aminoácidos, v) Efecto de la adsorción de aminoácidos en el 


transporte electrónico de una hoja de grafeno y, vi1) Conclusiones. 


Enfoque teórico/computacional 


Para el desarrollo de esta tesis se utilizaron algunos de los fundamentos teóricos 
aceptados en la actualidad para el modelado de materiales, como son la Teoría 
de la Densidad Funcional (DFT) junto con las Funciones de Green fuera del 
Equilibrio (NEGP). La combinación de estos formalismos teóricos, resulta ser 
adecuada para modelar sistemas que consisten en materiales conductores (o se- 
miconductores) entre dos electrodos sometidos a una diferencia de potencial. En 
el capítulo 2 de la presente tesis, se abordará con más detalle el enfoque teórico 
usado. El software implementado fue OpenMx, un código para el cálculo de la 
estructura electrónica de diversos materiales en la nanoescala, que ha dado exce- 
lentes resultados acordes con valores experimentales. Una de las ventajas de éste 
software es que permite realizar el modelado de sistemas para estudiar propie- 
dades de transporte electrónico, implementando el formalismo de DFT+NEGF. 
En el capítulo 3, se presentan los detalles computacionales y la descripción de 


los sistemas modelados. 


Secuencia de ejecución de la tesis 


Uno de los objetivos principales del trabajo, consistió en determinar los efectos 
de la adsorción de diez aminoácidos en las propiedades electrónicas del gra- 
feno: histidina (His), alanina (Ala), ácido aspártico (Asp), tirosina (Tyr), argini- 
na (Arg), glicina (Gly), ácido glutámico (Glu), fenilalanina (Phe), prolina (Pro) 
y lisina (Lys). Con el fin de obtener información completa, tanto del proceso 
de adsorción, así como también de sus efectos en las propiedades electrónicas 
del grafeno, el desarrollo de la tesis se llevó a cabo en dos etapas: 1) Resulta- 
dos y discusión de adsorción de aminoácidos y, 11) Efecto de la adsorción de 
aminoácidos en el transporte electrónico de una hoja de grafeno. En la primera 
etapa, se estudió la adsorción aminoácidos sobre una hoja de grafeno, se reportó 
valores para las energías y distancias de adsorción, configuraciones geométri- 
cas de equilibrio, densidades de carga y densidades de estado para cada sistema. 
En la segunda etapa, se calcularon las propiedades de transporte electrónico de 
una hoja de grafeno, antes y después de la adsorción de aminoácidos, donde 
se modeló un dispositivo de tres regiones: electrodo izquierdo (L), centro (C) y 
electrodo derecho (R). Los electrodos L/R corresponden a láminas de grafeno 
con 20 átomos de carbono cada uno. La región central de dispersión, contiene 
una hoja de grafeno con la molécula adsorbida (obtenida en la primera etapa). Se 
realizó el cálculo de la transmisión y de las curvas corriente-voltaje en una hoja 
de grafeno, antes y después de la adsorción de aminoácidos. En los capítulos 4 


y 3 se presentan las discusiones de los resultados para cada etapa. 


Datos significativos y hallazgos importantes 


Sobre el estudio de adsorción. Entre los principales hallazgos producto de la 
investigación realizada, se pueden mencionar los siguientes: 1) los aminoácidos 
se fisisorben sobre el grafeno a distancias inferiores a los 4 Á, donde, la naturale- 
za de sus cadenas laterales juega un papel preponderante en las configuraciones 
de equilibrio, así por ejemplo, los aminoácidos con cadenas laterales aromáticas, 
tienden a orientarse paralelamente a la hoja de grafeno. Este resultado es acorde 
a un trabajo reciente de autoensamblado de tirosina sobre grafeno realizado por 


Krebs y colaboradores; 11) durante los procesos de fisisorción el grafeno puede 


donar o aceptar electrones, en el caso de grafeno como aceptor, se modificó la 
naturaleza semi metálica del material. En todos los procesos de interacción se 
reportó la conformación de dipolos locales longitudinales y transversales entre 
la moléculas y los enlaces p del grafeno, esto concuerda con estudios experi- 
mentales realizados por Kumar Mallinen: y colaboradores; 111) las densidades de 
estados muestran la importancia de la inclusión de los grupos carboxilo y amino, 
ya que los estados cercanos a la energía de Fermi provienen en su mayoría de 


estos grupos funcionales. 


Sobre el estudio de transporte. En cuanto al estudio de los efectos de la fisi- 
sorción de aminoácidos en las propiedades electrónicas del grafeno, se pueden 
mencionar los siguientes aportes: 1) los perfiles de las curvas I-V son diferen- 
tes para cada aminoácido, lo que deja en evidencia la potencialidad del mate- 
rial en la construcción de dispositivos selectivos de sensado; 11) las curvas de 
transmisión se modifican para cada aminoácido dependiendo del potencial de 
polarización, así por ejemplo, aminoácidos como Phe y Asp, tienen una mayor 
transmisión para voltajes de polarización negativa, mientras que, la Lys permi- 
te mayor transmisión para voltajes positivos; 111) se observó una asimetría en 
los perfiles de las curvas de corriente-voltaje, lo que da lugar a otras posibles 
aplicaciones de aminoácidos en la construcción de nanodispositivos electrón1- 
cos basados en grafeno, como por ejemplo, desarrollo de diodos rectificadores 
moleculares; 1v) los aminoácidos Phe y Lys absorbidos sobre grafeno, presentan 
un comportamiento claro de diodo rectificador, en el primer caso, se inhibe el 
paso de corriente para voltajes de polarización positivos, mientras que se permite 


en valores negativos. En el segundo caso, se presenta el proceso inverso. 


Publicaciones y divulgación de resultados 


A lo largo de doctorado se publicaron cuatro artículos de revista, un capítulo de 
libro y, una parte de libro. Se participó en dos eventos científicos-tecnológicos 


con publicación. Se presentaron diez trabajos en eventos científicos. 


Premios y Distinciones 


Ganadora del concurso “Tesis en 3 minutos” de la Universidad Nacional del 
Litoral año 2017. 


Abstract 


Graphene is a bidimensional of carbon allotrope, with unique physical and chem- 
ical properties. The graphene is a flexible material, lightwe1ght, transparent, with 
high charge mobility and high surface area, among others. Since the experiment- 
ally obtainment of graphene in 2004 by Andre Geim and Konstantin Novoselov, 
the graphene has become in a promising material with applications in electronic, 


bioingner, biomedicine and nanotechnology. 


Particularly, in the biosensors building, the graphene is a material with numerous 
advantages: is sensitive, have a low electronic noise and 1s biocompatible. Sev- 
erals theoretical and experimental works have reported the ability of graphene 
to adsorb bacteria, cells and molecules, thus making 1t possible to study the 
eraphene properties before and after of adsorption process (and obtain inform- 
ation to functionalize and characterize the graphene-based biosensors). Works 


about of the DNA nucleobases have studied by several researchers. 


According to the interaction between graphene and biomolecules, the sensing 
mode can be: optical, magnetic, mechanical and electronic. In the literature, the 
field effect transistors based on graphene (GFETs) have gained great interest in 
recent years. The GFETS are electronic devices with three metallic lead (source, 
drain and gate) and one graphene chanel. When a voltage is applied in the drain 
(or gate), 1s possible evaluate of the electrical conductivity (or current) of the 


graphene before/after of the adsorption processes. 


In this work we have studied theoretically the effects in the electronics proper- 
ties of the graphene due to the adsorption process with ten amino acids. We have 
modeled a device with two graphene based-leads (source R and drain L) connec- 
ted to a graphene channel. According to the amino acid adsorbed on the channel, 
we can modeled the effect in the electric properties due to the charge transfer, the 


doping concentration and the formation of the local dipoles. We calculated the 


characteristic transmission and current-voltage curves (we considered positives 


and negatives polarization voltages: between -2 and 2 V). 


This work have divided in six chapters: 1) introduction, 11) theoretical concepts, 
111) description of the modeled systems, 1v) Results and discussion about the 
amino acids absorbed, v) Effects in the electronic transport of the a graphene 


sheet due to the amino acid adsorption and, vi) conclusions. 


Theoretical/Computational approach 


We used Density Functional Theory (DFT) and Non-equilibrium Green Func- 
tion (NEGEF) for the calculations in this work. The DFT+NEGF combination is 
the correct theoretical formalism for the conductor materials study (or semicon- 
ductors materials) between two electrodes. In the chapter 2 we approach with 


more details the theoretical framework. 


The implemented software was OpenMx. The methods and algorithms used in 
OpenMX has several advantages to investigate electronic structure properties 
of nanoscale materials, where the theoretical results are according to experi- 
mentally works. One of the advantages offered by OpenMx 1s the combination 
DFT+NEGEPF for the study of the electronic transport properties. In the chapter 3 


we presented the computational details about the modeled systems. 


Sequence of execution of this thesis 


One of the main objectives was study the effect on the graphene electronic trans- 
port properties due to adsorption of amino acids: histidine (His), alanine (Ala), 
aspartic acid (Asp), tyrosine (Tyr), arginine (Arg), glycine (Gly), glutamic acid 
(Glu), phenylalanine (Phe), proline (Pro) and lysine (Lys). With the aim to on- 
tained total information about the adsorption processes and graphene electronic 
transport properties, we have developed this work in two stages: 1) Results and 
discussion about to the amino acids adsorbed and, 11) Effects on the graphene 
electronic transport properties due to the adsorbed amino acids. In the first 
stage, we studied the amino acids adsorption on graphene and, we reported the 


adsorptions energies and distances, the equilibrium geometric configuration, the 


charge densities and the density of states for each system. In the second stage, 
we calculated the electronic transport properties in a layer of graphene before 
and after of the amino acid absorption. We modeled a device with three regions: 
left electrode (L), center (C) and right electrode (R). The L/R electrodes are 
graphene sheets with 20 carbon atoms each one. The central scattering region 
has a graphene sheet with adsorbed molecule (obtained in the first stage). We 
calculated (before and after of the amino acids absorbed), the transmission and 
the current-voltages curves. In the chapters 4 and 5 we reported the results for 


each stages. 


Significant data 


Adsorption study. The major research findings are: 1) the distance between 
amino acid and graphene sheet does not exceed the 4 Á, where, the side chain 
for each amino acid play a role important in the equilibrium configurations, e.g., 
the amino acids with aromatic side chains are oriented parallel to graphene sheet. 
This result 1s according to experimentally work by Krebs et al..; 11) the graphene 
can be an electron donor/acceptor. When, the Arg and Lys are adsorbed, changed 
the semi-metallic nature of the graphene. In all adsorption processes we re- 
ported the formation of longitudinal and transverse dipoles between molecules 
and graphene p-bonds, this result is according to experimentally work Kumar 
Mallinen1 et al. 111) the density of states show the important role of the inclusion 
of the carboxyl and amine groups, since the states near the Fermi energy are 


originated by this functional groups. 


Transport study. In the study of the effects of the amino acids physisorbed 
in the electronic transport in the graphene, we reported that: 1) the profiles of 
the curves I-V are differents for each adsorbed amino acid, this show the po- 
tentiality of graphene in the biosensors building; 11) the transmission changes 
for each amino acid and depends of polarization voltage , e.g. the Phe and Asp 
adsorbed amino acids have a major transmission when the polarization voltages 
are negatives, while that Lys has a major transmission when the polarization 
voltages are positives; 111) the curves profiles I-V are asymmetrics for the some 


molecules, this has application in the construction of graphene-based electronic 


nanodevices, e.g., in the rectifier diodes; 1v) the Phe and Lys adsorbed amino 
acids on graphene show the typical behaviour for the rectifier diode, Phe amino 
acid inhibits the currents when the polarization voltages are positive, while the 


Lys amino acid allows the current when the polarization voltages are negative. 


Publications and dissemination 


We have written four articles in journals, one book chapter and one part of the 
book. We have participated in two scientific events with publication and we 


presented ten works in scientist events. 
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Introducción 


En las últimas décadas, ha incrementado el interés por estudiar teórica y experimen- 
talmente, las interacciones generadas entre nanoestructuras basadas en grafeno y diferentes 
moléculas orgánicas, debido a las potenciales aplicaciones en el desarrollo de dispositivos de 
sensado en biomedicina, entre otros. Gracias a la alta sensibilidad eléctrica, propiedades físi- 
cas y químicas y, a la versatilidad para funcionalizar el material, el grafeno, es un candidato 


para la construcción de biosensores tanto a escala celular como molecular [|11|. 


1.1 Nociones básicas de Grafeno 


En el 2004 Andre Geim y Konstantin Novoselov lograron obtener experimentalmente gra- 
feno, una hoja de carbono de un átomo de espesor, estable a temperatura ambiente y con 
propiedades físicas y químicas únicas [111112]. 

Muchos investigadores habían estudiado por décadas los films ultra delgados y, desde 
sus experiencias, las monocapas eran imposibles de fabricar por diferentes razones: 1) ex- 
perimentalmente, evaporar películas metálicas de unos nanómetros de espesor daba como 
resultado formación de islas producto de la minimización de la energía en los sistemas y, a 
temperatura ambiente el material podría modificar sus propiedades físicas y químicas al reac- 
cionar con el aire —el oro por ejemplo, material inerte, se oxida parcialmente en el aire—, 


11) teóricamente el grafeno aislado podría ser termodinámicamente inestable. De acuerdo a 
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lo anterior, sintetizar por primera vez un cristal bidimensional, constituyó un gran descubri- 
miento para la época, Geim y Novoselov obtuvieron grafeno en principio, por exfoliación 
de grafito, usando cinta scotch para transferir y aislar una capa ultradelgada proveniente del 
material. Fabricaron contactos eléctricos con pintura de plata en los extremos de la capa y 
midieron algunas propiedades eléctricas. El grafeno exhibió alta conductividad, lo cual abrió 
el camino para la caracterización posterior de otras propiedades. En la actualidad hay dife- 
rentes técnicas para su obtención experimental y, aunque es más reactivo que el grafito, la 
interacción con el aire es débil, por lo que no influye en su red cristalina y no destruye su 
alta conductividad. En el año 2010 Geim y Novoselov obtuvieron el premio Nobel de Física 


por su descubrimiento. 


El grafeno es un alótropo de carbono con átomos dispuestos es una red hexagonal similar 
a un panal de abejas. La celda unidad contiene dos átomos ubicados en los sitios A y B, 
ver Fig. (1. 1[a). Los orbitales s y los orbitales degenerados Pp, y p, son los responsables 
de la hibridación sp? que da lugar a una estructura trigonal plana con formación de enlaces 
c, este tipo de enlace es el responsable de la fortaleza de la estructura de la red —con una 
separación entre átomos de carbono de 1.42 Á—. De acuerdo al principio de exclusión de 
Pauli, las bandas cd están llenas y conforman una banda de valencia profunda. El orbital p., 
perpendicular a la estructura planar, puede unirse mediante enlace covalente a otros átomos 
de carbono conformando enlaces T, cada uno con un electrón, por lo que la banda T está 


semillena (electrones deslocalizados). 


Este cristal bidimensional, es un semiconductor de gap cero (1.e., un semimetal) y sus 
propiedades electrónicas se originan principalmente por los electrones deslocalizados per- 
pendiculares a la estructura planar [113]. En los puntos de alta simetría K y K' de la primera 
zona de Brillouin del espacio recíproco (ver Fig. (1.1), la estructura de bandas presenta una 
relación de dispersión lineal cerca a la energía de Fermi, lo que conlleva a tener portadores 
de carga con masa efectiva cero “Fermiones de Dirac”. Entre las propiedades electrónicas 


del grafeno tenemos: 
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Figura 1.1: (a) Grafeno (b) Primera zona de Brillouin del grafeno (espacio recíproco) y (c) 
Bandas de energía teóricas para grafeno (se observa una dispersión lineal cerca a los puntos de 


alta simetría K y K”) 


e Bajo ruido electrónico y transporte balístico. En el grafeno, los portadores de carga 
a temperatura ambiente son débilmente dispersados, con una longitud de dispersión 


mayor a los 300 nm [114] 115]. 


e Excelente conductividad eléctrica (1738 S/m) [116]. De acuerdo a la conductividad 
medida en transistores de efecto de campo basados en grafeno, la conductividad incre- 
menta linealmente con el voltaje gate V, para las dos polaridades, donde la conducti- 
vidad nunca cae por debajo de un mínimo (valor correspondiente a la unidad cuántica 


de conductancia), incluso cuando las concentraciones de portadores de carga tienden 


a cero Fig. 1]. 


e Gran área superficial (2630 m' /g) y características ambipolares [117]. En el grafeno se 
puede controlar facilmente la densidad de portadores (electrones o huecos) mediante 


el uso de un campo eléctrico [18]. 
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e Alta movilidad de carga. Experimentalmente, se ha medido una movilidad de 15000 
em? Vs"? para una o tres capas de grafeno a temperatura ambiente, mientras que, 
para una hoja simple de grafeno a temperaturas cercanas a los O K la movilidad llega 
a los 230000 cm? V-*s7! (100 veces mayor que el silicio). 





pi Mi 


T=10K 


0.5— ca 











0 
V,(V) 
Figura 1.2: Conductividad, Adaptado de la Ref. 


La combinación de todas las propiedades electrónicas y su alta relación señal a ruido, le 
confieren al grafeno la capacidad de detectar cambios pequeños en las concentración de carga 
locales. Teórica y experimentalmente se ha mostrado que el grafeno interactúa con diferentes 
células, átomos o biomoléculas mediante fuerzas van der Waals y enlaces covalentes estables 
y, podría llegar a responder electrónicamente de manera diferente a cada compuesto, de tal 
forma que pueda usarse como un sensor altamente específico [20]. Adicionalmente, se han 
reportado cambios en el nivel de dopaje y en las propiedades electrónicas del grafeno en 
diferentes procesos de adsorción con diversas biomoléculas [9]. 

Estudiar y cuantificar los efectos electrónicos debidos a la interacción con diversos adsor- 
batos puede proporcionar herramientas para el desarrollo e implementación de dispositivos 
de detección basados en grafeno. Pero, ¿Cuáles son los fenómenos que subyacen a la trans- 
ferencia de carga entre el grafeno y las biomoléculas? ¿ Cómo funcionalizar el grafeno para 
hacerlo más selectivo y sensible a algunas biomoléculas? A continuación se expondrá algu- 


nas ideas básicas sobre la funcionalización del grafeno. 


1.2 Funcionalización: Grafeno como Biosensor 


1.2 Funcionalización: Grafeno como Biosensor 


Independientemente de la aplicación biomédica requerida, la funcionalización de cualquier 
material de sensado es importante para incrementar la selectividad y sensibilidad de los dis- 
positivos. En la literatura se consideran tres enfoques para mejorar la sensibilidad, especifici- 
dad y biocompatibilidad del grafeno a diversas biomoléculas o átomos: (1) funcionalización 
covalente, (2) funcionalización por dopaje y, (3) funcionalización no covalente [118|/211[22]. 

La funcionalización covalente es una técnica basada en la modificación de restos resl- 
duales sobre la superficie de grafeno (por ejemplo, el óxido de grafeno). Por otro lado, la 
funcionalización por dopaje, consiste principalmente en la sustitución de átomos de carbono 
de la red cristalina del grafeno. En la funcionalización (1) y (2), se altera inevitablemente la 
estructura electrónica original y las propiedades físicas del grafeno, modificando, por ejem- 
plo, los enlaces sp en sp3 , causando una disminución severa en la movilidad del portador 
[18]. Un tema importante sobre el dopaje consiste en controlar el tipo de dopaje, así como 
también la ubicación y concentración de dopantes. La lámina de grafeno puede doparse con 
electrones (dopaje tipo n) o agujeros (dopaje tipo p), métodos que han sido ampliamente 
estudiados [8]. 

La funcionalización no covalente se lleva a cabo por interacciones entre el grafeno y 
biomoléculas como ADN y péptidos, mediante fuerzas hidrofóbicas, van der Waals, elec- 
trostáticas y apilamiento T-7. Este tipo de funcionalización a diferencia de la covalente y 
por dopaje, resulta ser más conveniente porque las propiedades intrínsecas del grafeno orl- 
ginal se conservan. En este caso, se llevan a cabo principalmente procesos de fisisorción. 
En el cuadro. se presentan las características de dos procesos de adsorción que pueden 
presentarse entre grafeno y moléculas, átomos o impurezas. 

Diversos mecanismos basados en la funcionalización no covalente se han aplicado en el 
suministro de fármacos y genes, imágenes, ingeniería de tejidos y detección, ya que ofrece 
la posibilidad de adsorber de forma reversible moléculas en la superficie del grafeno sin 
alterar la red electrónica. Teóricamente, se han estudiado interacciones de diversas moléculas 
orgánicas con grafeno por sus potenciales aplicaciones en la construcción de dispositivos 
de sensado. En el cuadro. se sintetizan algunos resultados de las energías y distancias 
de adsorción de nucleótidos de ADN, aminoácidos y neurotransmisores. De acuerdo a los 
resultados: 1) hay una fisisorción fuerte entre nucleobases-grafeno y aminoácidos-grafeno, 


con distancias de adsorción superiores a los 2.8 Á, ii) el grafeno es un donor electrónico 
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débil y, 111) las interacciones no covalentes preservan la estructura hexagonal característica 


del grafeno. 


Cuadro 1.1: Comparación entre procesos de Adsorción física y química [8] [9] 


Proceso de | Energía de | Longitud de Modificación 
adsorción enlaces enlaces estructural 
Física 


10-100meV | >25A No 
Química | > 500meV 12.5 A Si 


Experimentalmente, se ha explorado la funcionalización por interacción con varias molécu- 








las policíclicas, derivados de pireno, ADN, aptámeros, fármacos aromáticos, colorantes, bio- 


moléculas y polímeros. Algunos trabajos se describen brevemente a continuación: 


e Las moléculas policíclicas tienen una afinidad con el plano basal del grafeno y el 
grafito, donde las interacciones de apilamiento 7-7 juegan un papel importante en la 
estabilidad de diferentes estructuras. Ghosh et al. funcionalizó el grafeno utilizan- 
do moléculas de carboxilatos (moléculas con una gran superficie plana y aromática y 
altamente soluble en agua), una fuerte interacción no covalente con el grafeno se ge- 
neró debido a la transferencia de carga. Esta funcionalización permitió estabilizar de 


grafeno (incluso de algunas capas) en soluciones acuosas. 


e Chen et al. confirmó la efectividad de la funcionalización no covalente del grafeno 
mediante la fabricación de la unión p-n de grafeno. En su trabajo, la característica 
tipo n o p para el grafeno bipolar, se moduló dopando una sola capa de grafeno con 
moléculas que donan o retiran electrones. Las moléculas 1,3 diaminonaftaleno (DAN) 
y 1-nitropireno (NP) fueron donadores de electrones (dopante de tipo n) y atrapadores 
de electrones (dopante de tipo n), respectivamente. Chen et al. reportó que los anillos 
aromáticos son responsables de la asociación de la molécula dopante con el grafeno, 


probablemente a través del apilamiento r-7. 
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Cuadro 1.2: Energías y distancias de adsorción para diferentes biomoléculas. Nucleobases: adenina 
(A), guanina (G), timina (IT), citosina (C), uracil (U); aminoácidos: histidina (His), alanina (Ala), ácido 
aspártico (Asp), tirosina (Tyr), valina (Val), arginina (Arg), asparginina (Asn) and cisteína (Cys); neu- 
rotrasmisores: dopamina (DA). La transferencia de carga del grafeno hacia las moléculas, en esta tabla 
se define como un valor positivo (negativo cuando las moléculas transfieren carga al grafeno). Los datos 
obtenidos en esta tesis se presentan y discuten con más detalle en el capítulo 4 —aquí sólo presentamos 
algunos datos calculados para His y Ala—.* distancia del átomo de carbono o hidrógeno más cercano a la 


superficie de grafeno. 


Ref Método Trans.Carga D“% (A) E“%s(eV) Biomolécula 
P5] MP2 — 3.50 0.94 
GGA-vdW 0.03 3.29 1.00 A 
DFT-vdW NN 3.5 0.63 
P5] MP2 — 3.50 1.07 
GGA-vdW 0.03 3.26 1.18 G 
DFT-vdW a 3.45 0.74 
MP2 == 3.50 0.80 
GGA-vdW 0.03 3.27 0.93 C 
DFT-vdW — 3.51 0.58 
MP2 — 3.50 0.74 U 
27) DFT-vdW — 3.49 0.54 
MP2 sun 3.5 0.83 
GGA-vdW 0.03 3.29 0.95 T 
P7] DFT-vdW — 3.53 0.60 
B3LYP-D3 = 3.36* 0.63 E 
DFT-D3 0.18 2.97 1.49 
B3LYP-D3 Es 2.63* 0.35 

Ala 
29]  DFT-D3 0.10 3.15 0.91 
[6] B3LYP = 3.21 3.79 e 
DFT-D3 0.17 4.00 1.17 dd 
[6] B3LYP an 3.07 3.89 Val 
[6] B3LYP a 3.14 4.21 de 
B3LYP-D3 a 3.05 * 1.66 i 
B3LYP-D3 > 2.65 * 0.51 Asn 
28] B3LYP-D3 a 3.45 * 0.44 Cys 
DFT-D2 0.25 3.18 0.74 

DA 
[31] GGA 0.01 2.84 0.46 


e Con el objetivo de detectar la actividad de fosfolipasa D, Si-Jia Liu et al. fun- 
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cionalizó el grafeno con fosfolípidos. Debido a la hidrofobicidad de la superficie de 
erafeno, las interacciones hidrofóbicas entre las colas de lípidos y el grafeno condu- 
jo un nanoensamblaje. En este trabajo se demostró que este nanoconjunto se puede 
adaptar a un novedoso biosensor de fluorescencia para ensayo de actividad de la fos- 
folipasa D mediante el uso de fosfoproteína marcada con fluoresceína lípido para la 


funcionalización. 


En síntesis, son muchos los fenómenos que subyacen en los procesos de funcionalización 
no covalente, como la interacción r-7 que estabiliza estructuras complejas, autoensambla- 
je molecular y la transferencia de carga hacia o desde el grafeno. Estudiando y analizando 
todas las interacciones generadas entre el grafeno y biomoléculas será posible desarrollar 
sensores de gran precisión con aplicaciones en control y diagnóstico de enfermedades para 
biomedicina. En la actualidad existen diversos métodos (eléctricos, electroquímicos, de fluo- 
rescencia, entre otros) para lograr el reconocimiento de biomoléculas de manera precisa, y 
en tiempo real. Dentro del enfoque de detección eléctrica, una posibilidad es la fabricación 
de transistores de efecto de campo basados en grafeno (GFETs) [33]. 


1.3 Transistores de Efecto de campo 


El dopaje de grafeno y la generación de dipolos eléctricos locales, son producto de los di- 
ferentes fenómenos que subyacen en los procesos de adsorción, tales como la transferencia 
de carga entre adsorbatos y sustratos, sin embargo, ¿es posible que la tranferencia de carga 
pueda modificar las propiedades electrónicas en el grafeno? De acuerdo a trabajos experi- 
mentales y teóricos, la movilidad electrónica del grafeno cerca de la superficie es extrema- 
damente sensible a los cambios de carga locales, por lo que varias propiedades electrónicas 
podrían alterarse en el proceso de adsorción —1ncluso cambios en el carácter semi metálico 
del grafeno— [34]. Uno de los dispositivos empleados para determinar los cambios en las 
propiedades electrónicas de grafeno, son los transistores de efecto de campo. 

Los transistores de efecto de campo basados en grafeno (GFET por sus siglas en inglés 
eraphene-based field-effect transistor) son dispositivos electrónicos con tres contactos metáli- 
cos (fuente “source”, drenaje “drain d” y compuerta “gate” ) y un canal de grafeno, los cuáles 
han sido estudiados para la detección de varias biomoléculas, ver Fig. 381. 


Al aplicar voltajes de polarización tanto en el drenaje d como en la puerta g (la fuente está 
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conectada a tierra y, el canal de grafeno y la compuerta están expuestos a un electrolito), 
es posible medir la respuesta del dispositivo, caracterizando la conductancia y las curvas de 
corriente antes y después de la interacción de la superficie del grafeno con biomoléculas 
[38]. Dependiendo de la biomolécula adsorbida, se pueden medir diferentes efectos en las 
propiedades eléctricas, debido a las diferencias en la transferencia de carga, la concentración 
de dopaje, el tipo de dopante y la generación de dipolos eléctricos locales. Además, en los 
GFET, las concentraciones de los portadores de carga se modifican por la presencia de un 


campo eléctrico externo, donde los electrones podrían reemplazarse por agujeros y viceversa 


[39140]. 






Puerta (Gate) 


Electrolito 


Drenaje 
EX (Drain) 


Figura 1.3: Dispositivo tipo GFET. 


En los dispositivos GFET el modo de detección es eléctrico, por lo que deben caracterl- 
zarse las curvas corriente vs voltaje (1-V) antes y después de que el grafeno haya adsorbido 
biomoléculas. Dos curvas típicas pueden llegar a analizarse: 1) corriente drenaje vs el voltaje 
gate (1 vs V,) y, 11) corriente vs tensión de polarización (6 voltaje bias V,). 

En las primeras, el V, tiene influencia en la conductancia a través de cambios en la mo- 
vilidad de carga [2]. Gráficamente, en una curva I vs V, (o curva de transferencia) el valor 
mínimo de conductancia corresponde al mínimo de transportadores de carga cerca al punto 
de Dirac (ver, Fig. (1.4)Ca,b)). A la izquierda del punto neutral, los transportadores son ma- 
yoritariamente huecos, mientras que a la derecha son electrones. El desplazamiento del nivel 
de Fermi se observa por desplazamientos de las curvas de transferencia —y consecuente- 


mente con el punto neutral— a lo largo del eje V,. Debido a la interacción con los diferentes 
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analitos, las curvas y sus respectivos puntos de Dirac pueden desplazarse horizontalmente, 


convirtiéndose así, en una herramienta de detección. 


(a) 
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Figura 1.4: (a) y (b) Curvas de tranferencia para GFET, (c) Curvas de corriente vs V; experl- 


mentales, para un dispositivo de detección bacteriana (Bacilus cereus) .. Adaptado de referencia 


y [3], respectivamente. 


En las segundas, se aplica una tensión de polarización V, entre dos terminales —diferencia 
de potencial entre la fuente y el drenaje V ps— [2//41]]. En esta configuración se omite el elec- 
trodo puerta o gate y, los cambios en la conductancia o resistencia del dispositivo, indican 
alteraciones físicas y químicas de las propiedades de la superficie (ver Fig. (1.4)(c)). 

En la actualidad, hay un avance en el desarrollo experimental y teórico de dispositivos 
GFET ya que, abarca la detección de señales a nivel celular e incluso molecular (biomolécu- 
las como la glucosa o nucleobases de ADN). Mohanty et al., fabricó y funcionalizó un 
dispositivo tipo GFET para la detección de la bacteria Bacilus cereus, en su trabajo se repor- 


tan cambios en las curvas I-V debido a la sensibilidad del dispositivo al exponerse a la bac- 
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1.3 Transistores de Efecto de campo 


teria (ver Fig. (1.4)(c)). Hess et al. [[4]] utilizando transistores de grafeno, registró potenciales 
de acción de células excitables, demostró citocompatibilidad en la interfaz neurona-electrodo 
y caracterizó la relación señal a ruido del dispositivo ver Fig. (1.5). Se reportó que el grafeno 
es un material prometedor para el diseño de nuevas prótesis neuronales, que tendrían como 
función no sólo el registro de actividad eléctrica sino la estimulación [42]. Por otro lado, a 
escalas aún menores, pruebas experimentales de sensores químicos de nanotubos de carbono, 
para moléculas pequeñas como NO», NH, Oz y Ha han demostrado la alta sensibilidad en la 


estructura electrónica del grafeno cuando adsorbe diferentes moléculas en su superficie [43]. 


contactos metálicos 
célula 





a a a 


"A a 
e a e 


Grafeno 


Figura 1.5: Dispositivo para registrar señales a nivel celular, adaptado de referencia 


En el cuadro (1.3), se listan algunas propiedades de los transistores de efecto de campo 
con diferentes materiales propuestos para aplicaciones en bioelectrónica. El grafeno exhibe 
una biocompatibilidad comparable con el diamante y el AIGaN/GaN, aunque cuenta con ma- 
yor movilidad de carga — 33 veces más que el diamante y 3 veces más que el AIGaN/GaN— , 
capacitancia interfacial'| y flexibilidad [110]. Además, los GFET tienen notables ventajas so- 
bre dispositivos equivalentes basados en nano tubos de carbono; son fáciles de fabricar, y la 
naturaleza ambipolar del grafeno tiene una capacidad inherente para responder tanto a los 


adsorbatos cargados negativa como positivamente [45]. 


!Se puede entender como la propiedad de almacenar carga eléctrica cuando hay una diferencia de potencial 


entre dos superficies, por unidad de área. Se tiene que 7 = tra + E (donde e es la constante dieléctrica 





del material aislante), la capacitancia depende del espesor dieléctrico, cuando d —> O se tiene la capacitancia 
interfacial [44]. 
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1. INTRODUCCIÓN 


Cuadro 1.3: Propiedades de algunos transistores de efecto de campo construidos con diferentes 
materiales (Seifert M. et al.) [10]. Se toma como referencia la biocompatibilidad y la flexibilidad 


del silicio (denotado con Ref.), los símbolos + y - indican mayor y menor, respectivamente. 









Movilidad | Capacitancia 
Material de carga y | interfacial | Biocompatibilidad | Flexibilidad 
(cm/Vs) 
Silicio 
Diamante 


Teóricamente, Lee at al. modeló un dispositivo tipo GFET para estudiar los efectos 
de la adsorción de nucleótidos de ADN (Adenina (A), citosina (C), guanina (G) y timina (T)) 
en la conductancia de nanocintas de grafeno. El modelado se llevó a cabo a través de cálculos 
de primeros principios y funciones de Green fuera del equilibrio. Los resultados sugieren que 
la detección de A, C, G y T puede ser monitoreada por los cambios en la conductancia de la 
nanocinta de grafeno luego de adsorber cada nucleótido. En esta misma dirección, cambios 
en las curvas I-V de nanocintas de grafeno debido a la adsorción de nucleótidos fueron 
reportados por Song et al. [46]. 

Los dispositivos GFETs son muy versátiles, por lo que aún es necesario estudiar diferen- 
tes mecanismos de optimización, como por ejemplo, estudiar la maximización de los efectos 
de las biomoléculas en los canales de conducción (adaptando la geometría del dispositivo), o 
explorar el control de la selectividad y sensibilidad del dispositivo mediante la manipulación 
de la conductividad del grafeno. 

En este trabajo proponemos el modelado teórico de un dispositivo tipo GFET's para eva- 
luar su eventual sensibilidad a diferentes aminoácidos y, valorar su potencial aplicación como 


secuenciador de proteínas. 


1.4 Objetivos de la Tesis 


Las proteínas de todos los seres vivos son indispensables ya que, cumplen una gran variedad 


de tareas en el funcionamiento de cada organismo y además, todos los procesos biológicos 
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1.4 Objetivos de la Tesis 


dependen de la presencia o la actividad de este tipo de moléculas —desde funciones estructu- 
rales, tales como el colágeno que integra las fibras altamente resistentes en tejidos de sostén, 
hasta funciones de transporte en la sangre como la hemoglobina—. 

Las proteínas se encuentran conformadas por largas cadenas de aminoácidos y cualquier 
modificación en la secuencia de los mismos, afecta directamente la correcta función de la 
proteína. Conocer la secuencia de aminoácidos es de relevancia ya que: 1) brinda informa- 
ción sobre su función, estructura e historia —mecanismos de acción—, 11) posibilita elaborar 
nuevas proteínas con nuevas funciones (sólo modificando el orden de los aminoácidos), 111) 
contribuye en el conocimiento de las reglas que gobiernan el plegamiento de las cadenas 
polipeptídicas, 1v) juega un rol importante en el estudio de enfermedades genéticas, una al- 
teración en la secuencia de tan sólo un aminoácido tiene consecuencias nefastas, como por 
ejemplo, la anemia falciforme (producción anormal de glóbulos rojos) o, la osteogenia im- 
perfecta (también llamada huesos de cristal) donde ocurre el cambio del aminoácido glicina 
por una cisteína en una posición específica de la proteína|!]. 

De acuerdo a lo anterior, resulta útil pensar en la fabricación de dispositivos que permitan 
conocer rápidamente la secuencia de aminoácidos en una proteína. Aprovechando las pro- 
piedades del grafeno en el desarrollo de diversos mecanismos de detección eléctrica a escala 
molecular, en este trabajo se estudió teóricamente los efectos que produce la fisisorción de 
aminoácidos en las propiedades electrónicas del grafeno, para evaluar su eventual aplicación 
en un secuenciador de proteínas. Para ello, se modeló mediante funciones de Green fuera 
del equilibrio (NEGF por sus siglas en inglés Non-Equilibrium Green's Function), un dis- 
positivo tipo GFET con dos electrodos de grafeno y, se calcularon curvas de corriente vs 
voltaje bias (1-V;,) para determinar si el dispositivo respondía eléctricamente diferente con 


cada aminoácido (1.e., si es selectivo a cada aminoácido). 


Determinar la secuencia de un aminoácido en las proteínas permite por un lado caracterizarlas y detectar 
anomalías moleculares responsables de diferentes enfermedades y, por el otro, controlar la calidad de proteínas 
producidas en ingeniería genética. En 1945 W.H. Stein y S. Moore realizaron el primer análisis completo de 
aminoácidos de una proteína la b-lactoglobulina, utilizando métodos químicos y microbiológicos. Estas técni- 
cas, y otras equivalentes, fueron utilizadas por diferentes grupos de investigadores. Sin embargo, determinar la 
secuencia completa de una proteína era una tarea que demandaba varios años. CH Li, por ejemplo, demoró diez 
años para establecer la de la hormona de crecimiento humana. A partir de 1950 el método de degradaciones de 
Edman ha sido utilizado para secuencia cadenas largas, donde los grupos terminales se etiquetan y se separan 


del péptido sin afectar los enlaces. Una de las limitaciones de éste método, es que no es 100 % confiable. 
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1. INTRODUCCIÓN 


1.4.1 Objetivos Generales 


e Modelizar el comportamiento electrónico de diversas estructuras basadas en grafeno 


mediante teoría de la funcional densidad y funciones de Green fuera del equilibrio. 


e Simular la adsorción de moléculas biológicas en grafeno para aplicaciones en nano 
dispositivos 


e Estudiar del transporte cuántico en nanosistemas para aplicaciones en nanobiosenso- 


res. 


1.4.2 Objetivos Específicos 


e Estudiar la fisisorción de aminoácidos sobre una superficie de grafeno mediante cálcu- 
los de primeros principios, para determinar las configuraciones de equilibrio, energías 
y distancias de adsorción, densidades de estado y transferencia de carga entre substra- 


tos y adsorbatos. 


e Modelar un dispositivo tipo GFET implementando funciones de Green fuera del equi- 
librio, para estudiar los efectos que produce la fisisorción de aminoácidos en las pro- 


piedades de transporte electrónico del grafeno. 


e Evaluar la posibilidad de su usar grafeno en la construcción de un secuenciador de 
proteínas mediante el análisis de la curvas de transmisión y de corriente vs voltaje de 


polarización. 


El modelado se llevó a cabo considerando la adsorción de diez aminoácidos: histidina 
(His), alanina (Ala), ácido aspártico (Asp), tirosina (Tyr), arginina (Arg), ácido glutámico 
(Glu), glicina (Gly), fenilalanina (Phe), prolina (Pro) y Lisina (Lys). En el capítulo tres se 
expondrá en detalle la estructura interna de los aminoácidos comtemplados en el estudio. 

En el anexo B, se presentan los trabajos publicados en el marco del desarrollo de ésta 


tesis. 


'En el anexo A, se presenta otra aplicación de la adsorción de moléculas sobre grafeno, mediante: Auto- 
ensamblaje molecular sobre grafeno. Esta estrategia de funcionalización, permite estudiar nuevos materiales 


para el desarrollo de sensores y dispositivos de almacenamiento de energía. 
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Conceptos teóricos 


El modelado computacional del dispositivo tipo GFET, se desarrolló dentro del enfoque 
de la teoría de la funcional densidad en combinación con funciones de Green fuera del equi- 
librio. En este capítulo se describirán los fundamentos y conceptos teóricos sobre los cuales 


basamos los cálculos ab initio en el presente trabajo. 


2.1 Teoría de la Funcional densidad 


La teoría de la funcional densidad (DFT de sus siglas en inglés “Density funcional theory”), 
es un método de la mecánica cuántica, para resolver sistemas de muchos cuerpos que incluye 
la correlación electrónica (interacciones entre los electrones) en los cálculos de la estructura 
electrónica del estado fundamental de moléculas y sólidos. 

El estudio de muchos cuerpos, conlleva a considerar N electrones en presencia de M 
núcleos, examinando a su vez, cómo se organizan las partículas de tal forma que perma- 
nezcan estables en el tiempo. Desde la mecánica cuántica significa resolver en principio, la 


ecuación de Schródinger dependiente del tiempo. 


OE, E5),0) 


Ol = Hrv(ter), (En), t) (2.1) 
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2. CONCEPTOS TEÓRICOS 


La autofunción 4((€7), 1€),t) es dependiente del tiempo t y de las coordenadas espa- 
ciales y de espín de todas las partículas, tanto electrones (£;), como núcleos [€;). Try es el 


operador hamiltoniano total: 




















M N M N 
» hi hi 1 | —lel(Z,lel) 
Hr = VS — NE A MIE 
A) td) E ko 

N ON M M 

l l [A 1 1 |(Zalel Zgle!) 

e A JUAel 
2293) ÁTEO ==] 22D ATEO Ro — Rol 


donde, M denota cantidad de núcleos, siendo R. La y Ma la posición, carga y masa nuclear 
respectivamente. N hace referencia a la cantidad de electrones, r; y m, a la posición y masa 
del 1-ésimo electrón. 

Cada uno de los términos de la ec. corresponde a un tipo de interacción entre las 


partículas que componen el sistema: 


e f., es el operador de energía cinética de los núcleos 





Ñ a 
Elen ES 3 a ve) (2.3) 





No 2 
É.., = Y DA 
2, E j (2.4) 
e V,, es el operador de energía potencial electrón-núcleo 


=lel(Zalel) 


o Ra 


(2.5) 


— 2 2 == 


a=1 1¿=1 








e V.. es el operador de energía potencial electrón-electrón 





a la 1 [Ele 
V=3 22 2 0 o e” 


2.1 Teoría de la Funcional densidad 


AN 


e V,,, es el operador de potencial núcleo-núcleo 





M M 
, 1 1 |(Zalel)(Zalel) 


El operador hamiltoniano total l p[!]se escribe de forma más compacta: 


MAN AN 


Hr 3% an eN Evo un e on Vo. + Ven (2.8) 


La ec. es imposible de resolver por métodos análiticos y numéricos, por lo que es 
necesario: 1) considerar que el hamiltoniano total es independiente del tiempo, de tal forma 
que, la función de onda se puede expresar como el producto de dos funciones, una que de- 
penda de las coordenadas espaciales (electrónicas y nucleares) y otra, de las coordenadas 
temporales, siendo Yr (LE), 1€5),8) = Vr(LEr+, 1€7))x(t). Esto permite separar variables 
y obtener una ecuación independiente del tiempo y, 11) recurrir a la aproximación de Born- 
Oppenheimer para núcleos fijos?] y, solucionar la ecuación de Schródinger para electrones y 
núcleos por separado. Además por considerar los núcleos fijos, el término Vs, del hamilto- 
niano electrónico tendrá un valor constante. Lo anterior permite orientar el problema y resol- 
ver una ecuación de Schródinger independiente del tiempo, para núcleos fijos, donde la au- 
tofunción electrónica dependerá de 3N coordenadas y de las interacciones electrón-electrón, 
y no de todas las coordenadas espaciales, de espín y temporales de núcleos y electrones de 
la ec. (2.1). En síntesis, los estudios teóricos de las propiedades electrónicas estáticas inde- 


pendientes del tiempo, se reducen a resolver la ecuación exclusivamente para electrones: 


'Una primera aproximación para resolver la ecuación de Schródinger dependiente del tiempo, consiste en 
no incluir interacción entre espines e interacciones espín-órbita, por lo que, el hamiltoniano total H7 no las 


considera. 
“Teniendo en cuenta que el movimiento del núcleo es menor en comparación con los electrones, y que la 


razón entre la masa del protón y el electrón da cuenta de un núcleo muy masivo, 1:1836, Born y Oppenheimer 
en 1927 proponen separar el movimiento de electrones y núcleos. Mediante cálculos semiclásicos, comparan las 
escalas de energías de 1) un electrón que tiene una transición entre dos estados sucesivos, 11) energía rotacional 
del núcleo y, 111) energía vibracional del núcleo. Se encuentra que la energía del electrón es mucho mayor que 


la energía vibracional y rotacional del núcleo. 
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2. CONCEPTOS TEÓRICOS 





Donde V.... de la ec. es el potencial electrostático generado por los núcleos de la 
ec. (2.5), por notación en este documento se asume Vo == Vos. 

La ecuación de Schródinger independiente del tiempo de la ec. (2.9), es válida para N 
electrones que interactúan con M núcleos fijos y puede ser resuelta a través de diferentes 
aproximaciones. Cuando el enfoque del problema consiste en hallar inicialmente la mejor 
función de onda 4(r;,73,-+- ,Ty) para describir el sistema de partículas, se trabaja desde 
la aproximación de Hartree-Fock. Si el tratamiento se centra en la densidad electrónica p(r”) 
como variable elemental, se está situado en el formalismo de la DFT. 

Muchos observables experimentales dependen de la densidad electrónica, lo que permi- 
te pensar que el cálculo de 4(7;,73,--- ,ry) arroja más información de la requerida en la 
práctica. Uno de los primeros trabajos en esta dirección, fue desarrollado por L.H. Thomas 
y E. Fermi (1927), quienes calcularon todas las propiedades de interés, partiendo de la den- 
sidad electrónica [47!. 

Basados en esta idea Pierre Hohenberg y Walter Kohn demuestran que la densidad electróni- 
ca p(r) del estado fundamental determina unívocamente un potencial externo V..,, suficiente 
para construir todo el operador Hamiltoniano de la ec. (2.9), y finalmente obtener todas las 
propiedades del estado fundamental [48]. En síntesis, p(r) da cuenta de todas las propieda- 
des del estado fundamental, o en otras palabras, las propiedades del estado fundamental son 
funcionales de p(7). 

En 1964 Hohenberg y Khon formularon las ideas esenciales de la DFT mediante la de- 


mostración de dos teoremas por reductio ad absurdum (reducción al absurdo) [49]. 


e Teorema 1. El potencial externo V..,. es (a menos de una constante) un único funcional 
de p(r). Dado que V....; fija el H, es posible ver que el estado fundamental de todas las 


partículas son funcionales únicas de p(7). 


Corolario: Todas las propiedades del sistema están completamente determinadas por 


p(r). 
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2.1 Teoría de la Funcional densidad 


e Teorema 2. Se puede definir una funcional universal de energía E[p| válida para cual- 


quier sistema de N electrones. 


De acuerdo al teorema 1 de Hohenberg y Khon, la energía y el número total de electrones 
N, que son propiedades del estado fundamental, son a su vez funcionales de la densidad 


electrónica 


E =Elp] (2.10) 


IN = Jona (2.11) 
Además cualquier densidad de prueba ¿(7) que obedezca la relación de la ec. (2.11), 


cumple 


Elp] > Elp (2.12) 


lo que implica a su vez que existe un principio variacional para el funcional de energía 
verdadero. 

Kohn y Sham (KS) en 1963 reformularon el problema, al dejar en evidencia que no es 
necesario calcular la función de onda completa de N cuerpos v(r;,73,-+- ,Ty) para deter- 
minar las propiedades del estado fundamental de un sistema de partículas, basta con hallar 
la densidad electrónica p(r) [48]. En 1998 Kohn gana el premio Nobel en química por el 
desarrollo de la Teoría de la Funcional Densidad. 

Para definir las bases de la DFT, Kohn y Lu Sham escribieron la energía electrónica total 
como funcional de p(7), y resolvieron una ecuación tipo Schódinger para partículas simples. 


Dos de los señalamientos más importantes en esta teoría son: 


1. La energía electrónica total de la ec. (2.9) es una propiedad fundamental, por tanto, es 


posible expresarla como un funcional de la densidad electrónica, así: 


El =T19+ $ 05 Voto + Vat] d+ Ela] 2.13) 


donde, T'[p] corresponde a la energía cinética, f p(P) Ves (7)dr energía debida al po- 
tencial externo, 5 f Ve(F)p(F)dF interacciones electrón-electrón y É.,. energía de inter- 


cambio y correlación. 
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2. Con el fin de simplificar los cálculos y reducir la solución de la ecuación de Schródin- 
ger a N ecuaciones de una partícula simple, Kohn y Sham introducen el concepto de 
cuasi electrones, que son ferminones no interactuantes que se mueven en un poten- 
cial externo V.,,, construido a partir de la densidad total del sistema real, es decir, la 


densidad electrónica del sistema real, es igual a la de los cuasi electrones [47]. 


En este orden de ideas, al aplicar el principio variacional para la energía total como 
funcional de la densidad, llevar a cabo un estudio de propiedades electrónicas en el marco de 


la DFT, implica resolver las ecuaciones de Kohn y Sham (KS) para cuasi electrones: 


(2.14) 





siendo V. y ,, el potencial efectivo de KS deducido a partir de la comparación de los resultados 
de aplicar el principio variacional para el sistema real y el de cuasi electrones" | Matemática- 


mente la expresión para el potencial efectivo es de la forma: 


Very (r) = Venir) + Volr) + Veclr) (2.15) 
donde 


e V..,:(r) es el operador de interacción núcleo-electrón. 


' Aplicando el principio variacional E[p] > El[p], con la restricción N = f p(r)dF a través de multiplica- 
dores de Lagrange, se tiene para el sistema real: 


5le=a(frmar-x)|=o 


ST dEzc 
+ Vent?) + Vo(F = 
59 +Ver Vol += 
Y para los cuasi electrones 
ST 
A 
o prkr) =p 


Donde el multiplicador de Lagrange y es por construcción el potencial químico. Al igualar las expresiones del 


sistema real y cuasielectrones, se define Voy (7) = Veas (7) + Volr) + Sas E En = | 
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0 Va(r) es el operador de interacción electrón-electrón. Puede obtenerse desde la ecua- 


ción de Poisson dadas las apropiadas condiciones de frontera 


VVal(r) = —4re*p(r) (2.16) 


e V.. (7) término de correlación e intercambio. 


E 
AS de (2.17) 


mu 


Exec = Exc + (1 —T) (2.18) 


En DFT se incluye el término del potencial de correlación electrónica e intercambio, el 
primero hace referencia a la interacción coulombiana electrón-electrón, mientras que 
el segundo corresponde a la definición dada desde la aproximación de Hartree-Fock. 
El potencial de intercambio y correlación V..., y la energía de intercambio y corre- 
lación £,,., integran además en forma aproximada la diferencia de energía cinética 
entre el sistema de N electrones interactuantes 7”, con el de N cuas1 electrones no- 


interactuantes 1”. 


El estado fundamental para este sistema se encuentra dado por la ocupación de las N 
partículas simples en los orbitales energéticamente más bajos, y la densidad electrónica total 


está dada por: 


po() => wi (51 (2.19) 


Por construcción, esta densidad es igual a la densidad electrónica del sistema real de 
interés, sin embargo, los orbitales de partículas simples 4, y las energías e, no son en princi- 
pio definidos para los electrones del sistema real. A pesar de ello los resultados de cálculos 
realizados por DFT demuestran que constituye una muy buena aproximación. 

La energía total de los cuasi electrones es la suma de las energías para partículas simples 
e, de los orbitales ocupados, sin embargo, la energía total del sistema de electrones es 1gual 


d: 
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2. CONCEPTOS TEÓRICOS 


N 


Elp| = y (il — 


¡=1 
Mediante la ec. y la ec. (2.15) la energía electrónica total del sistema se puede 


escribir como: 


apo) + [ at0) [Vaar) +3Vo(0) a(7) + End] (220 


2m 


N 


Elp] = DE = ] p(r) ¿Val A Vr) d(r) + Exclp! (2.21) 


Es importante resaltar, que el hamiltoniano de las ecuaciones de KS queda determinado 
a partir de pp, tal y como se indica en la ec. (2.14), donde todos los términos son funcionales 
de la densidad de carga. Sin embargo, el hamiltoniano incluye variables de intercambio y 
correlación V... y E... que no tienen forma exacta como funcionales de p(r) . Esto conlleva 
a seleccionar alguna aproximación desde el formalismo de la DFT, para enunciar comple- 
tamente la energía electrónica total. Cabe señalar, que si se conociera la forma exacta de 
Ver) y Exc, se tendría la solución verdadera para la energía y propiedades fundamentales 
del sistema en estudio. 

Entre las diferentes aproximaciones para el potencial de intercambio y correlación se 
tienen la aproximación densidad local LDA y gradiente generalizado GGA. Existen otro tipo 
de aproximaciones meta-GGA, híbridos, y DFT que incluye interacciones tipo van der Waals 
(DFT+vdW), entre muchas otras. 


2.1.1 DFT-D3: GGA semiempírico 


El origen de la interacción de van der Waals entre dos fragmentos que no están unidos quím1- 
camente, es producto del acoplamiento del campo eléctrico generado por las fluctuaciones en 
la densidad electrónica entre fragmantos. Este efecto en la correlación dinámica, no es captu- 
rado por los funcionales locales y semilocales habituales, tales como LDA y GGA, y no está 
relacionado con el intercambio, de modo que el £,, por sí solo tampoco ayuda. A grandes 
distancias, las interacciones de van der Waals se modelan como interacciones dipolo-dipolo 
clásica, que decae como C,/ RR”. El potencial de Lennard-Jones es ampliamente utilizado 
para representar las interacciones no enlazantes en estos contextos [150]. 

Las interacciones débiles no covalentes responsables de la complejidad de las macro- 


moléculas o interacciones intramoleculares, corresponden a los enlaces de hidrógeno (o 
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puentes de hidrógeno) y fuerzas de van der Waals (vdW). Estas últimas incluyen interac- 
ciones dipolo-dipolo, dipolo-dipolo inducido y fuerzas de dispersión de London. Interaccio- 
nes de este tipo entre átomos y moléculas juegan un papel importante en muchos sistemas 
químicos y, junto con la interacción electrostática y de intercambio-repulsivo controlan la 
estabilidad, por ejemplo, de estructuras como ADN, proteínas, y formación de agregados, 
entre otras. 

Las vascos de intercambio y correlación LDA o GGA-PBE, no describen com- 
pletamente vdW [51] [52]. Por tal motivo, se llevaron a cabo cálculos teniendo en cuenta 
interacciones de leo alcance, usando la aproximación semiempírica de S. Grimme. DFT- 
D3 [54], que implementa el programa OpenMX. 

La idea general de la DFI-D3 propuesta por Grimme en el año 2010, consiste en reem- 
plazar los efectos de correlación electrónica no local, de largo y medio alcance del gradiente 
convencional, por términos que dependen de un coeficiente de dispersión por par de electro- 
nes C,, R7”. La base fundamental está en las series de Becke desarrolladas en 1997 [55]. 


La energía total será de la forma: 


Eprr-D = Exs-DrT + Edisp (2.22) 


siendo Exs-prr la energía obtenida del autoconsistente de las ecuaciones de Kohn y Sham 


para DET, y Eqisp es la corrección de dispersión empírica. 


Edisp = Y 2, > E famnol Ras) (2.23) 


er: n=6,8,10,. 

La suma sobre todos los pares de átomos del sistema. Donde, 02? es la media (isotro- 
pica) del coeficiente de dispersión de orden n (ordenes n = 6,8,10,...) para cada par de 
electrones AB [.J nmémol"?]. Ras es la distancia interanuclear. s,, es un factor de escala 
que depende del funcional densidad. fdamp es una función de amortiguación para pequeños 
valores de R. Los valores empíricos de Cé y los radios de van der Waals se listan en el trabajo 


desarrollado por S. Grimme [541]. 
2.1.2 El método de campo autoconsistente (SCF) 


Las ecuaciones de Kohn y Sham, ec. y ec. (2.19), se resuelven mediante el método de 
campo auto consistente (SCP), en la Fig. se presenta un esquema del método de SCF, 
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para la solución de las ecuaciones de KS. Tenemos: 


Selección de la Estructura 


Densidad de carga 
inicial 





Q(F) 


¿qe A 
op Y + Visio) Y (A =P) 






Auto consistente? 
Converge? 





Variables de interés 


Figura 2.1: Ciclo autoconsistente (SCHP). 


e Se parte de una densidad de carga inicial py(r). 


e Se calcula el V.,,, ya que es funcional de po(7) (adoptando una aproximación de co- 


rrelación e intercambio). 


e Conocido V. ;,, el hamiltoniano queda definido para resolver las ecuaciones de KS$. 
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2.1 Teoría de la Funcional densidad 


e La solución de KS arroja resultados para autofunciones Y y autovalores e, que sirven 


para hallar una nueva densidad electrónica p; (7) 


e Con la nueva p(7) se re-calcula el potencial efectivo V. y y y se repite el proceso, hasta 
que las densidades de entrada y de salida sean iguales ( su diferencia igual o menor a 


un cierto valor de tolerancia) considerando el cálculo convergido. 


e Cuando se obtenga el valor para p(7) convergido, se calculan todas las propiedades 
que dependen de esta magnitud (DOS, bandas de energía, etc) y que son de interés 


para el sistema en estudio. 


2.1.3 Solución de Ecuaciones de KS para partícula simple 


Los métodos DFT para cálculo de estructura electrónica, se basan en la expansión en alguna 


base de las eigenfunciones 4, de la ecuación de KS, así 


N; 
r) = 3 Cijój(F) (2.24) 


donde /,(r) son las funciones base, C,; los coeficientes de expansión y N;, el número de 
bases p,(7). La densidad electrónica, también se puede escribir en términos de esta expansión 


de la siguiente forma: 


N N, 


=> Y CAC (Por (r) (2.25) 


¿=1 7, k=1 
Una vez se determinan las bases ¿;, se definen los elementos de la matriz de overlap $ y 


los elementos de la matriz del hamiltoniano H: 


Six = (Ojl0x) (2.26) 


h? 


Diagonalizando la ecuación secular: 
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se obtienen las energías e, de KS y, los orbitales de las ecuaciones ec. (2.14) y ec. (2.24). 
Usando ahora ec. se re-calcula Vo y V..., y se sigue el procedimiento SCF. 

Existen diversos métodos DFT que se diferencian unos de otros, por el tipo de repre- 
sentación que realizan de los potenciales y las funciones base. El desarrollo de los métodos 
tiene como objetivo minimizar el costo computacional y de programación, manteniendo la 
precisión en los cálculos. Entre los métodos más conocidos se tienen el LAPW (Linearized 
Augmented Plane Wave), el APW (Augmented) y el de pseudopotenciales. 


En este trabajo se adopta el método de pseudopotenciales. 


2.1.3.1 Pseudopotenciales 


Los estados electrónicos de un átomo se clasifican en 1) estados de core, que son los electro- 
nes que se encuentran altamente localizados cerca a los núcleos atómicos y no participan en 
los enlaces químicos, 11) estados de valencia, son los electrones de las capas externas sem1- 
llenas de los átomos, responsables de los enlaces químicos y, 111) estados de semi-core, son 
localizados y polarizables pero no están directamente relacionados a los enlaces químicos 
[SO]. 

Los electrones del core, interactúan fuertemente con el núcleo y no responden con efi- 
ciencia a los electrones de valencia, por lo que podrían llegar a considerarse fijos. Por otro 
lado, los electrones de valencia son sensibles al entorno atómico, lo que conlleva a que la 
densidad de carga varíe para el electrón que está aislado en su estado fundamental. Además 
los electrones de valencia se ven afectados por los electrones de core mediante un potencial 
externo. 

Con el fin de agilizar los tiempos de cálculo y reducir el tratamiento matemático a sólo 
los electrones de valencia, el efecto del core sobre los electrones de valencia, se represen- 
ta mediante un pseudopotencial (aproximación de core congelado), que no es más que un 
potencial de interacción que actúa sólo sobre los electrones de valencia. De esta manera, se 
eliminan estados centrales y oscilaciones de ortogonalización en las funciones de onda de 
valencia, ver Fig. (2.2)[5]. 

Los pseudopotenciales son una buena aproximación de la función de onda a partir del 
radio de corte r. (medido desde el núcleo). Para que describan correctamente el interior de 
la esfera de radio r,, se construye un pseudopotencial que reproduzca una pseudofunción de 
onda que no tenga nodos en el interior de la esfera, pero que conserve la norma, es decir, que 


a pesar de que la densidad de carga verdadera y la pseudodensidad de carga no sean 1guales 
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en la región del core, la carga integrada sobre dicha región sea la misma en los dos casos, se 


tienen pseudopotenciales ab initio conservadores de la norma. 





yPS 


Figura 2.2: Esquema de pseudopotenciales y pseudofunciones de onda p*” (líneas pun- 


teadas). Las funciones reales y las pseudofunciones coinciden a partir del radio de corte r,. VP* 


y pP* carecen de nodos en la región del core [5]. 


2.2 Funciones de Green Fuera del Equilibrio 


En esta sección se presenta una introducción de los conceptos teóricos relacionados al es- 
tudio del transporte electrónico a escala molecular, donde las leyes físicas que gobiernan la 
conducción eléctrica se describen haciendo uso de la mecánica cuántica. 

Teóricamente existen diferentes enfoques para abordar problemas de transporte cuántico: 
1) Métodos semiempíricos, los cuáles no son autoconsistentes y son basados en parámetros 
tieth-binding para expresar el Hamiltoniano y el sistema molecular aislado. 11) Método de 
superceldas, las cuáles son basadas en la solución de las ecuaciones de Khon-Sham con con- 
diciones periódicas de contorno. Los estados de dispersión o de scattering son determinados 
mediante técnicas recursivas (este método no describe sistemas con diferentes electrodos ni 
los sistemas bajo un potencial bias externo). 111) El método jeliumm de Lippman-Schwinger, 
donde los electrodos son descritos en términos del modelo jellium y las ecuaciones de Khon- 
Sham son resueltas autoconsistentemente para la estructura abierta. La densidad de carga se 


describe desde las regiones de dispersión del dispositivo. Este método no tiene en cuenta los 
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estados ligados que existen dentro del dispositivo (sólo usa los estados de dispersión para 
construir la densidad de carga y el potencial). 1v) Las funciones de Green fuera del equili- 
brio (NEGF por sus siglas en inglés non-equilibrium Green's function) [56]. Este esquema 
ofrece muchas ventajas ya que la densidad de carga es determinada desde la teoría propuesta 
por Keldysh [57], además permite estudiar esquemas típicos de dispositivos con dos elec- 
trodos conectados a una región central (tratados por el mismo marco teórico) y, se pueden 
introducir efectos de múltiples cuerpos, e.g., electrón-fonón o electrón-electrón, en forma de 
autoenergías sin desviarse sensiblemente del marco teórico [58]. 

En la presente tesis, el estudio se dirige hacia la descripción de la estructura electrónica 
y las propiedades de transporte de uniones moleculares con electrodos de grafeno, dichos 
sistemas no son finitos ni periódicos, lo cuál hace complejo aplicar sólo la teoría de la DFT 
—además, el sistema está fuera del equilibrio cuando se aplica un voltaje bias—. De acuerdo 
a lo anterior, la teoría más conveniente para abordar el problema de transporte, es la comb1- 


nación DFT+NEGPF. 


A continuación se explicará brevemente el modelo teórico. 


2.2.1 NEGF para un sistema con electrodos 


El estudio de transporte, se enfoca en el modelado de un canal de grafeno acoplado a dos 
electrodos del mismo material, en este caso, se asume que los electrones de la hoja de grafeno 
se transmiten a través del sistema sin sufrir dispersión inelástica (es decir, el transporte se 
da en el régimen balístico). Una vez el grafeno ha adsorbido una molécula o es dopado, los 
electrones que viajan a través del sistema experimentan una dispersión, incrementando la 
probabilidad de que algunos electrones no puedan ser transmitidos [56]. 

En este orden de ideas, el formalismo teórico se basa en el transporte elástico, es decir, 
en ausencia de excitaciones vibracionales. En el régimen de acoplamiento fuerte —que es 
el caso a tratar en este trabajo—, los orbitales de la molécula se hibridizan con los de los 
contactos, y los electrones pueden viajar de forma coherente de un electrodo a otro. Por lo 
tanto, los estados propios de la molécula aislada ya no presentan una base adecuada para la 
descripción del transporte, siendo la corriente dominada por los procesos de tunel elásticos, 


donde no existen los modos de vibración. Debe aclararse que, en la implementación de NEGF 
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en este trabajo, no incorpora los efectos de decoherencia por eventos de scattering inelásticos 
(pérdida de la memoria de fase por interacción con fonones, o con los eletrodos) 

La coherencia del movimiento de los electrones, que no se observa en los sistemas ma- 
Ccroscópicos, es notoria en nanoestructuras y, surgen varios fenómenos interesantes, como 
la cuantificación de la conductancia, las características inusuales de tensión de corriente no 
lineal y el efecto de bloqueo de Coulomb. 

Abordar el estudio del transporte electrónico en la nanoescala desde el método de las 
NEGFs tiene múltiples ventajas, ya que permite: a) analizar desde el mismo marco teórico 
todas las partes que constituyen un dispositivo de dos electrodos conectados a una región 
central; b) determinar la estructura electrónica de la región central auto consistentemente, 
combinando cálculos ab initio por DFT o Hartree-Fock; c) introducir los efectos de múltiples 
cuerpos (como por ejemplo, electrón-fonón, electrón-electrón) mediante auto energías. A 
continuación, se realiza una breve descripción del fundamento teórico implementado por 
Ozaki et al 160]. que se enfoca en el transporte coherente y que es fundamento de los 
resultados presentados en el capítulo cinco de ésta tesis. 

El problema central consiste en describir el transporte electrónico de dispositivos de dos 
electrodos conectados a una región central, que se someten a una diferencia de potencial entre 
los electrodos. Este problema se considera fuera del equilibrio y, puede abordarse perturbat1- 
vamente partiendo de la solución para el problema sin potencial aplicado. En la Fig. se 
representa el esquema de un sistema general compuesto de tres celdas: izquierda (L, ), central 
(Co) y derecha (R;). Para construir la función de onda ye de todo el sistema (zona central y 


electrodos) se usa una base de pseudo potenciales atómicos (PAOSs) (4: 


1 
Y (r) = 75 e ria Pia (Y — Ti — Ru) (2.29) 


Donde d denota el espín (UP o DOWN), yu es el índice de banda, R,, = /¿b + /¿c es el 
vector de red, k= kb + k,c es el vector de onda de Bloch (con b y e vectores base de la 
(k) 


TU 


red en el plano yz, y b y € sus recíprocos), los números c son los vectores a determinar 
auto consistentemente, los vectores 7, dan las posiciones atómicas, el índice «: barre todos 
los orbitales utilizados para describir cada átomo. De acuerdo al método descrito, se adopta 
periodicidad en el plano yz, quedando x libre de esa restricción. Considerando la variación 


1Si se requiere estudiar los efectos de decoherencia en NEGF, un modelo a estudiar, puede ser el desarro- 
llado por D/Amato-Pastawski [59] 
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Figura 2.3: Configuración de un sistema típico para ser tratado por el método de NEGF. La 
periodicidad del sistema se da en el plano yz, donde la región central C es una región extendida 


que contiene Co, Lo y Ro. Los electrodos R y L son infinitos a lo largo del eje x. 


total de la energía con DFT de un sistema expresado por las funciones de KS con respecto a 


los coeficientes c, se obtiene la ecuación matricial KS. 


HN) = SO (2.30) 


gp 
(k) 
gu, 0 


k 
Donde cS o es un vector columna con los coeficientes [c 


», la matriz del Hamiltoniano 
A Lo. y la de solapamiento (que se debe a la no-ortogonalidad de los PAOs), 5%, se definen 


como: 


a e y SS (2.31) 


n 


SA = Sy At Ra (2.32) 


n 
Donde hosiajóR, Y Siajór, SOn los elementos de matriz del Hamiltoniano y la matriz 
de solapamiento entre dos orbitales, d,(r — 7) y Pjalr — 7; — R,,). Ahora se considera 


una región central extendida C compuesta por Co, Lo y Ry como se muestra en la figura. 
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Adicionalmente se imponen dos condiciones: 1) los PAOs, (/, de la región C( se solapan con 
los de las regiones Ly y Ro, pero no con los de las regiones L; y R;; y 11) los PAOSs de la 
región L;(R;,) se superponen sólo con los de las regiones vecinas más cercanas, L;+¡(R;+1). 


Estas condiciones permiten escribir la matriz del Hamiltoniano como: 


a , 
k k 
k Mad, sigo (K) 
NS Ho ss Be (2.33) 


Donde los elementos A e H Do y H he: . Son matrices cuadradas del tamaño de la can- 
tidad de orbitales utilizados (Nz,, No,N,.,) para la descripción de las regiones, L;¡, C y Ry, 
respectivamente. Los elementos fuera de la diagonal conectan las distintas regiones. La ma- 
triz de solapamiento tiene una estructura análoga (aunque cabe remarcar que S'“% no de- 
pende del espín, d). Adaptando una de las definiciones usuales de la función de Green, 
G(Z) = (Z— H)* (donde Z = E + ie es una variable compleja con su parte real, E, 
representa una energía), al modelo aquí planteado escribimos EL ) = (Z5W -— qe) 


y, desagregando el Hamiltoniano en su elemento central, H a, y el resto, escribimos: 


ODE TDS =D (2.34) 


Donde la influencia del electrodo izquierdo y derecho es tratada como una perturbación 


de la región central (C) mediante auto energías 2 Z ) y iZ ): 


k k k k k k 
AZ) = (286%, - Hi CS: — Hito) (2.35) 


k k k k k k 
210%) A (A Mi IC Ol > NE Hna) (2.36) 


O 


Donde CU(Z ) y az ) son las funciones de Green de los electrodos derecho e iz- 
quierdo. Por las condiciones 1) y 11) mencionadas arriba, las matrices del Hamiltoniano y de 


solapamiento para las superficies izquierda y derecha también son trigonales: 
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A: HH: 0 

(ki) Ha Ho Haz 

H% = Ho Ha (2.37) 
a | 
911 912 0 

(K) 931 992 933 

su) — So Se (2.38) 
0 


donde s= L, R. El subíndice 1 que aparece en los elementos H;,, (o S;¿) se refiere a la 
celda L, (o R;). Todas las celdas L;(1 = 1,2,3,...) son idénticas entre sí, al igual que las 
de R, entre sí (aunque L y R pueden ser diferentes entre ellos) y son calculadas (densidad 
electrónica, hamiltoniano y matriz densidad) por separado de la región central, por un cálculo 
estándar con condiciones periódicas en 3D, es decir, sin tener en cuenta la región central. Esta 
suposición es aceptable si la región central se diseña de modo tal que cerca de lo frontera 
con los electrodos, la diferencia sea mínima, por tal motivo se debe recurrir a la región C 


extendida. 


2.2.1.1 Matriz densidad en equilibrio 


Si el Hamiltoniano y la matriz de solapamiento asociadas a la ec. son dependientes de 
k, esto implica que la función espectral de Lehmann de la función de la región central es un 
número complejo. Entonces la matriz densidad oo donde uno de sus orbitales asociados se 
encuentra en la celda central y otro en la celda ubicada en R,,, se puede expresar por medio de 
funciones de Green retardada y avanzada, a, (E+10%) ff Ep) y au, (E-10%)f(E—-p), 


respectivamente, según: 


e dl _i 
E A O (2.39) 
Cc yBZ 


con 


pe => ] dEGUA(E 440%) f(E — p) (2.40) 
e 


siendo V. el volumen de la celda unidad, f 7 la integración sobre la primera zona de 


Brillouin, 0* un infinitesimal positivo, f es la función de distribución de Fermi-Dirac y u el 
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potencial químico. La integración en el espacio k se desarrolla numéricamente por medio de 
una red tipo Monkhorst Pack [61]. 
Si el Hamiltoniano y las matrices de solapamiento son independientes de k, se puede 


simplificar la matriz densidad de acuerdo a la expresión conocida: 


dEG/,(E+i05fAE =p) (2.41) 





La matriz densidad posibilita el cálculo de la densidad de carga de la región central (C), 


que a su vez, se descompone en tres contribuciones: 


E SI EA A (2.42) 
s=L,R 


donde, s = L/R,ne“(r), y“ (r) y n*s (7) son las contribuciones a la densidad de carga 
provenientes de las funciones base cuyos orbitales se encuentran ubicadas en la región cen- 
tral (C), electrodo y centro (sC) y electrodos (ss), respectivamente. El factor 2 del segundo 
término de la ec. proviene de considerar la contribución 1“*(r). Explícitamente cada 


contribución se define como: 


7 Sá (a => y o 10,10,R,, Dio (Y e T;)O55(r o 15 o R,,) (2.43) 


n 10,18 


ner => y Pot ¡QRr,¡BR,, ¡Oia Y — Ti > (Ra se a) ds (r AA R,,:) (2.44) 


nn! 10,1B 


A = A Y rr ar ¡Dial — 7 — (Ra +8)]ójo[r — 7, -(R ta) (45) 


nn! 10,18 

donde a es el vector de red de la celda unitaria para las regiones Ly y Kio a lo largo del 
eje x. El desplazamiento —a(+a) denota que las función base está en L, (RR, ) en la dirección 
del eje x. Las densidades Do Br,, Y Poor, ¡gr,, SON calculadas por DFT convencional 
con condiciones periódicas 3D (sin que sea necesario aplicar NEGP), siendo aplicadas a 
los electrodos por separado antes de realizar el cálculo de las funciones Green de la cadena 
infinita representada en la Fig. (2.3). Consistentemente, el potencial de Hartree en la frontera 
entre la región central y los electrodos se considera igual a la de los electrodos, es decir, que 
se toma el potencial de frontera C' — L(R); igual a la de la frontera L(R)¡ — L(R)a. 
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2.2.1.2 Matriz densidad fuera del equilibrio 
Basados en la teoría de NEGF desarrollado por Schwinger y Keldysh [[57]], la matriz densidad 


fuera del equilibrio de la región central, se evalúa como: 


PoR. = Pos E APO) (2.46) 
Se adiciona a la matriz densidad en equilibrio pa un término de corrección APO Ri) * 
Se define: 
1 
Apora => | deApre E (2.47) 
Ve JBz 
siendo 
1 
ApW = e 'N dECCUE + iO (EGUE —iAf(E) (2.48) 
Tis 
con 
TP (E) = 420 (E+ io —-EN (E —ie)) (2.49) 
y 
Af(E) = FE — Us) Aa) (2.50) 


El potencial químico más bajo (sea el yr O LL) se usa para calcular la matriz densidad 
en equilibrio. Las ec. y ec. se calculan teniendo en cuenta que sí; = Ry sz = L 
Si UL < PR y, $1 = e Ri un < < HL. 

En el contexto de las NEGF se deben tener en cuenta dos suposiciones para resolver 
adecuadamente la ec. (2.46): 1) la ocupación de las funciones de onda incidentes desde la 1z- 
quierda/derecha obedecen a la función de distribución de Fermi Dirac de dichos electrodos, 
incluso cuando ya están en la región central; 11) en la región central los estados con energía 
por debajo del potencial químico menor, están en equilibrio debido a otros obstáculos físi- 
cos tales como interacciones electrón-fonón y electrón-electrón los cuáles no se consideran 
explícitamente en la implementación utilizada. La última suposición permite que los elec- 
trones que se encuentren por debajo del menor potencial químico puedan ocupar estados 
altamente localizados por medio del primer término de la ec. (2.46). Sólo los estados en 
régimen de energía entre los dos potenciales químicos son tratados como en no equilibrio en 


los cuáles la contribución de las funciones de onda incidentes desde el electrodo con mayor 
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potencial químico son tomados en cuenta para formar el término de correción dado por la 
ec. (2.47). 

Finalmente, en la ec. (2.48), el integrando no es analítico ya que es función de Z y Z*, 
por lo que no se puede aplicar el método de integración de contorno utilizado en la matriz 
en equilibrio. En este caso se utiliza un esquema de cuadratura rectangular simple en el eje 
real con una pequeña parte imaginaria e. Dado que el integrando contiene una diferencia 
entre dos funciones de distribución de Fermi-Dirac, el rango de energía para la integración 


se reduce ampliamente a un pequeño intervalo. 


2.2.1.3 DFT+NEGF 


Resolver sistemas de transporte cuántico en el marco de referencia de DFT combinado con 


NEGF puede resumirse en el siguiente bucle iterativo auto consistente [162]: 


n"(r) > DFT > Hxs[n(r)] > NEGF > n*(r) (2.51) 


El cálculo parte de una densidad de estados electrónico inicial n'” (r) que se emplea para 
definir el hamiltoniano de Kohn-Sham Hxs|n(r)| mediante cálculos DFT. Posteriormente la 
inversa de la matriz permite definir la función de Green retardada cuya integración sobre la 
energía determina la matriz de densidad a través de la fórmula basada en NEGF. La matriz 
de elementos n““(r) = (r|p|r) será la nueva densidad electrónica para empezar una nueva 
interacción. Este proceso se repite hasta alcanzar el criterio de convergencia Á, que converge 
cuando se cumple la condición: p"* — p*'” < A. Una vez los resultados convergen es posible 
obtener la transmisión y la corriente a través del dispositivo de dos electrodos. 

La transmisión es evaluada por la fórmula de Landauer para la región central C 


conectada a dos electrodos: 


1 
T(E) == / dae?Tw(p) (2.52) 
Ve az 
Donde 
TO(Ey) =Tr [TE (EGUAE +00 (EJGUA(E — ie) (2.53) 


En el esquema de Landauer-Búttiker, la corriente / a través de un conductor se expresa 


en términos de la transmisión, según: 
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dE - dET(DTAE) (2.54) 





La teoría de Landauer-Búttiker es hoy día un formalismo ampliamente utilizado para 
el cálculo de las propiedades de transporte electrónico de sistemas en la nanoescala. Este 
formalismo es apropiado para el estudio de los sistemas con electrodos propuestos en el 


presente trabajo. 


36 


Descripción de los sistemas 
modelados 


En este capítulo se presenta por un lado, algunas nociones básicas sobre la estructu- 
ra general de los aminoácidos estudiados en esta tesis y, por el otro, la descripción de los 
parámetros computacionales que se implementaron en las dos etapas del presente trabajo — 
1) estudio de la adsorción de diez aminoácidos sobre grafeno y, 11) estudio de los efectos de 


la adsorción de aminoácidos en las propiedades electrónicas del grafeno—. 


3.1 Aminoácidos 


Los aminoácidos son moléculas que conforman a las proteínas. En el genoma humano exis- 
ten 20 aminoácidos que son los que componen las proteínas y por lo tanto son llamados 
proteinogenicos (ó alfa-aminoácido). La estructura general de un alfa-aminoácido compren- 
de un carbono central (alfa) unido a un grupo carboxilo (-COOH), un grupo amino (-NH»), 
un hidrógeno y una cadena lateral (o cadena radical R)|'| ver Fig. (3.1). La cadena lateral 
es específica para cada aminoácido, determinando la identidad y las propiedades de cada 


molécula. 


lalfa-aminoácido, también hace referencia a que los grupos funcionales carboxilo y amino están unidos al 


mismo carbono. 
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(A) Carboxilo 


Figura 3.1: Esquema de la estructura general de Aminoácidos. Grupos funcionales y cadena 


lateral. 


El estudio de adsorción de moléculas sobre grafeno se llevó a cabo para diez aminoácl- 
dos: Histidina (His), Alanina (Ala), Ácido Aspártico (Asp), Tirosina (Tyr), Ácido Glutámico 
(Glu), Glicina (Gly), Fenilalanina (Phe), Arginina (Arg), Prolina (Pro) y Lisina (Lys). De 
acuerdo a la cadena lateral que los caracteriza, las moléculas consideradas se agrupan de la 
siguiente manera: 1) aminoácidos ácidos o con carga negativa a pH fisiológico, como en el 
caso de Asp y Glu, éstas moléculas tienen un segundo grupo carboxilo en su cadena secunda- 
ria; 11) aminoácidos aromáticos como Phe y Tyr, que poseen cadenas laterales aromáticas de 
benceno y fenol, respectivamente; 111) aminoácidos básicos (cadenas laterales muy polares) 
como Lys, Arg e His; iv) aminoácidos alifáticos que tienen una cadena lateral hidrocarbona- 
da (o en otras palabras, constituidos por carbono e hidrógeno cuyo carácter no es aromático), 
como la Gly y Ala. La prolina también tiene una cadena lateral de naturaleza alifática, pe- 
ro difiere de los demás en que su cadena lateral está unida tanto al carbono alfa como al 
nitrógeno del grupo amino. En la Fig. y Fig. se presentan los aminoácidos con su 
cadena lateral característica. 

En este trabajo, el modelado de los procesos de adsorción se desarrolló en fase gas —1.e., 
no se contempló la ionización de los grupos funcionales debido a cambios de pH, puesto que, 
el aminoácido no se encuentra en solución 166] —, por lo que, se mantienen los grupos 
funcionales -COOH y -NH» en todo los aminoácidos. 


La discusión detallada sobre la importancia de incluir los grupos funcionales en este es- 
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3.1 Aminoácidos 





Figura 3.2: Aminoácidos considerados en este trabajo. (Derecha) Los átomos de la cadena lateral 
se resaltaron con color verde en esta figura, con el fin de facilitar al lector el reconocimiento de 


las cadenas laterales carcteríticas de cada molécula. 


tudio, se llevó a cabo en las referencias —Atrabajos publicados por los autores de 
la presente tesis—. Los cálculos teóricos muestran diferencias significativas para energías de 
adsorción y distancias de adsorción, entre el aminoácido con y sin grupos carboxilo y amino. 
Considerar el estudio de la biomolécula con -COOH y -NH» es relevante, ya que, tal como 
se mostró para la His, su presencia desplaza los estados que provienen del anillo más cerca a 
la energía de Fermi del grafeno, lo que favorece la eventual funcionalización y aplicación del 
grafeno en la construcción de biosensores —se resportan diferencias al comparar las curvas 


E-V de la His con y sin grupos carboxilo y amino [69] —. 
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Figura 3.3: Aminoácidos considerados en este trabajo.(Derecha) Los átomos de la cadena lateral 
se resaltaron con color verde en esta figura, con el fin de facilitar al lector el reconocimiento de 


las cadenas laterales carcteríticas de cada molécula. 


3.2 Detalles computacionales 


Para las dos etapas propuestas en este trabajo —estudio de adsorción y de transporte electrónico— 
, los cálculos ab initio se llevaron a cabo en el marco de la DFT, mediante la aproximaciones 
de intercambio y correlación GGA-PBE, con correcciones de dispersión de largo alcance 
en el esquema DFT-D3 |'| Se usaron pseudopotenciales conservadores de la norma del soft- 


ware OpenMx y se implementó pseudo-orbitales atómicos (PAO) como funciones base para 


¡Como se discutió en el capítulo 2 y, de acuerdo a trabajos teóricos de adsorción de moléculas, el método 
DFT-D3 mejora la exactitud de los cálculos, demostrando ser consistente con resultados experimentales para 
sólidos y superficies [70)[711. 
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expandir la ecuación de Khon-Sham. Un PAO se especifica como H5.0-s1, C5.0-s2p2dl, 
06.0-s2p2d1, N6.0-s2p2d, donde H, C, O y N corresponden a los símbolos atómicos y los 
valores 5.0, 6.0 al radio de cut-off de la función PAO en unidades de a.u. Además, la nota- 
ción s2p2dl significa que hay dos estados s, dos estados p y un estado d en la función radial 


60) 1721. 


A continuación se presentan las superceldas para el estudio de adsorción y de transporte. 


3.2.1 Caracterización de la supercelda para adsorción 


El sistema substrato-adsorbato constituido por un aminoácido y una hoja de grafeno, no 
corresponde propiamente a un cristal, sin embargo, todas las herramientas y aproximaciones 
que se usan para el estudio de la estructura electrónica de sólidos, pueden implementarse sin 
problema, ya que, todos los cálculos que se llevan a cabo desde la DFT arrojan información 
para los electrones de la celda primitiva en la primera zona de Brillouin (1BZ). En este 
caso, para el sistema en consideración, sólo basta definir una celda unidad que contenga en 
su interior todos los átomos de la hoja de grafeno y de la biomolécula, en otras palabras 
una supercelda. Tenemos así, una “estructura cristalina” compuesta por traslaciones de la 
supercelda, aunque los cálculos obtenidos corresponden a la celda unidad. Para caracterizar 
correctamente la supercelda y evitar la interacción entre moléculas de celdas vecinas, fue 
necesario tener en cuenta los siguientes cálculos previos: 


1) Para hoja de grafeno 


e Convergencia de puntos k. De acuerdo a cálculos de trabajos anteriores [[68l, para 
una hoja de grafeno de 60 átomos de carbono, con dimensiones 12,78 Áx 12.2975 


Áx 25 Á, la grilla óptima de puntos kes de 5x5 x 1. 


e Energía de corte (cut-off). En la Fig. se graficó la energía total en función a la 


energía de corte. La energía total converge para valores de 180 y. 


e Distancia entre capas de grafeno desacopladas. En la Fig. (3.5) se graficó la energía 
total del grafeno puro en función de la separación entre capas (parámetro c), las capas 


se encuentran desacopladas a partir de los 6 Á . 


11) Para los aminoácidos 
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Figura 3.4: Energía total vs. Energía cut-off. Grafeno. 
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Figura 3.5: Energía total vs. separación entre capas de grafeno (parámetro c). El desacople entre 


capas se puede observar a partir de los 6 Á. 


e Es conveniente establecer el valor del lado de la celda cúbica para el cual no hay inter- 
acción entre moléculas de celdas vecinas, para ello se elaboró una gráfica de energía 
total de la molécula en función del tamaño de la celda —+en este caso particular, pre- 


sentamos los resultados para His—. De acuerdo con la Fig. (3.6), la energía total de la 
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molécula converge a partir de una celda de tamaño 10Áx 10 Ax 10Á. 
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Figura 3.6: Energía total vs. Lado de celda cúbica. Histidina. 





Figura 3.7: Supercelda histidina sobre grafeno. Los átomos de oxígeno, carbono, hidrógeno y 


nitrógeno, se representan con color rojo, café, rosado y gris, respectivamente. 


Siguiendo los resultados anteriores, el sistema total sustrato-adsorbato se estudió consi- 
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derando una celda ortorrómbica simple de lados 12,78 Áx 12.2975 Ax 25 Á, cuya base es 
una supercelda de grafeno de 60 átomos de carbono. Se consideró una energía cut-off de 180 
Ry y una grilla k de 5 x 5 x 1, siendo el refinamiento de la grilla para el calculo de la DOS 
de 50 x 50 x 1 puntos k. La distancia de separación entre aminoácidos y grafeno es h, ver 
Fig. (3.7). El criterio de convergencia para la relajación fue de 0.02 eV/Á (fuerza máxima 


entre átomos) para todos los aminoácidos. 


3.2.1.1 Metodología para determinar la energía y distancia de adsorción 


La energía de interacción (Ej) y la distancia de adsorción (d.,¿,) entre substrato y adsorbato, 


se calculó para cada molécula, de acuerdo al siguiente procedimiento: 


l. Primera Relajación: Se relajó la molécula —sin substrato— para eliminar tensiones 


internas 


2. Segunda Relajación: La molécula relajada, se ubica a una distancia de 4.5 Á de una 
hoja de grafeno|*| Se efectúo una segunda relajación de todo el sistema —átomos de 
la molécula y de la hoja de grafeno—. Se exploró diferentes orientaciones para la 
molécula con el fin de determinar la configuración de mínima energía —la Ey de- 
pende de la disposición atómica—. 

A manera de ejemplo, en la Fig. (3.8), Fig. se presentan las configuraciones 
estudiadas para ALA y ASP, respectivamente. Á partir de todas las orientaciones estu- 


diadas para cada aminoácido, se selecciona la configuración de menor energía. 


3. Energía de interacción en función de la distancia entre substrato-adsorbato: Se ubicó 
la molécula —cuya disposición geométrica dió como resultado una minimización de 
energía— sobre la hoja de grafeno a diferentes distancias, efectuando una relajación 
completa del sistema. Se obtienen dos datos para cada distancia considerada: 1) la 
distancia media final /h de la molécula (distancia promedio del centro geométrico de la 

La estructura 3D de las moléculas se encuentran disponibles en la base de datos Protein Data Bank PDB 


http://www.rcsb.org/pdb/home/home.do [73] 
“De acuerdo a estudios teóricos y experimentales, los aminoácidos y las nucleobases del ADN, interactuán 


con el grafeno a distancias menores a 5 Á[67] [74]. 
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Figura 3.8: Alanina sobre grafeno (a) grupo -COOH enfrentado a la hoja de grafeno (b) grupo 
-NH> enfrentado a la hoja de grafeno (c) grupos -COOH y -NH» enfrentados a la hoja de grafeno. 














Lo, (a) LL, (b) Ly (0 


Figura 3.9: Ácido Aspártico sobre grafeno (a) grupos -COOH y cadena lateral enfrentados a la 
hoja de grafeno (b) grupo -COOH enfrentado a la hoja de grafeno (c) grupos -COOH y -NH» 


enfrentados a la hoja de grafeno. 


molécula en la dirección Z con respecto al centro geométrico del grafeno en 21!]y, 11) 
la energía total E**(h) para cada configuración. Con los datos anteriores se elaboró 


una gráfica de energía de interacción E”"*(h) en función de la altura media final h. 


La energía de interacción se define como: 


ER NN pt NN pref (3.1) 


Medir la distancia entre la molécula y la hoja de grafeno, tiene algunas variantes. En el proceso de ad- 


sorción, la hoja de grafeno se deforma, de tal forma que los átomos de carbono que se encuentran justamente 


debajo la molécula tienden a alejarse o acercarse (dependiendo de las fuerzas de atracción o repulsión). Los 


átomos del grafeno que se encuentran lejos de la molécula, es decir, los que se encuentran en los extremos de 


la base, no se mueven significativamente en el eje z. Por tal motivo la distancia h se define como la distancia 


entre el centro geométrico en Z de la molécula y el centro geométrico del grafeno en z. El tiempo de cálculo 


por punto es de 4 a 6 días (cálculos en paralelo en cluster Coyote y Aquiles). 
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donde E"*% es una energía de referencia, que equivale a la energía total para cuando 
la distancia de separación entre substrato-adsorbato es de h = 12,5 Á (el grafeno y 
la molécula se encuentran lo suficientemente alejados y no se presenta interacción). 
La energía de adsorción Eds es el mínimo global de la curva —en la Fig. 
se presenta un esquema cualitativo de la curva de interacción— |*| que resulta de la 
ec. (3.1). 


Distancia entre 
substrato-adsorbato 


Energía interacción 





! 

y / 
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Aa 


Ni 
Distancia de Adsorción 


Figura 3.10: Esquema cualitativo de la curva de energía de interacción vs. distancia entre 


substrato-adsorbato. 


En el capítulo 4 se discuten las energías y distancias de adsorción calculadas para cada 


sistema. 


3.2.2 Caracterización de la supercelda para cálculos de Transporte 


Para el tratamiento de propiedades de transporte electrónico, en este trabajo se estudió un 
dispositivo tipo GFET como el presentado en la Fig. (3.11). El dispositivo se compone de una 
celda central Cy que contiene los átomos de una hoja de grafeno y un aminoácido adsorbido 

| La adsorción se modelizó realizando una relajación de todos los átomos que componen el sistema para 
cada altura considerada. En la literatura se proponen otros métodos que no relajan completamente el sistema 
—o relajan parcialmente—. Se llevó a cabo una comparación entre métodos, lo que permitió ratificar la impor- 


tancia de una relajación total para encontrar el mínimo global de las curvas, con el fin de no recaer en mínimos 


locales. En la referencia [67|, se presenta con más detalle la discusión. 
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—adoptando las configuraciones de equilibrio del cálculo de adsorción que se llevó a cabo 
en la primera etapa del trabajo—. La celda extendida C, la cual incluye las celdas de los 
electrodos Ly y Ry, tiene dimensiones de 21.30 Áx 12.29 Áx 25.00 Á. Las celdas Lo.1.2.. y 
Ro. 1.2... tienen dimensiones de 4.26 Áx 12.29 Áx 25.00 Á con 20 átomos de carbono cada 


uno. 


Vo 














Figura 3.11: Supercelda para cálculos de transporte. (a) esquema de supercelda, (b) vista lateral, 
(c) vista superior. La región C extendida comprende las regiones: Cy (determinada en la primera 


etapa del trabajo), Lo y Ro. 


El método de NEGF, permite calcular las corrientes a través de los cálculos SCF a todos 
los voltajes de polarización, sin embargo, la selección adecuada de los voltajes, depende 
de la naturaleza de los sistemas en estudio. El modelado que se llevó a cabo en esta tesis, 


comprendió un rango de voltajes bias entre -2 V y 2 V, tal y como ha sido reportado en otros 


trabajos teóricos 


!Si bien, este trabajo constituye un modelado teórico —lo que permite en principio considerar voltajes bias 
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mayores—, no es adecuado hacer un modelado con valores bias mayores a 2 V, ya que, experimentalmente, las 


nanocintas de grafeno presentan ruptira eléctrica alrededor de los 1.8 V 
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Resultados y discusión de 
Adsorción de aminoácidos 


En este capítulo se presentan los resultados del modelado de la adsorción de diez aminoáci- 
dos sobre grafeno. Se discuten las configuraciones más estables, las densidades de estado y 
la transferencia de carga entre substratos y absorbatos. Para facilitar la lectura de los resulta- 
dos, en los próximos capítulos se nota His/Grap, Ala/Grap, Asp/Grap, Tyr/Grap, Arg/Grap, 
Glu/Grap, Gly/Grap, Phe/Grap, Pro/Grap y Lys/Grap, para hacer referencia a la adsorción 
sobre grafeno de un aminoácido en particular (1.e., histidina adsorbida sobre grafeno, se nota 


como His/Grap). 


4.1 Resultados de Configuración más estable 


Las configuraciones de equilibrio para los diez aminoácidos fisisorbidos se presentan en las 
Fig. y Fig. (4.2). 

De acuerdo a las geometrías más estables, se observan diferentes disposiciones de los 
grupos funcionales (-COOH y -NH») en los proceso de adsorción, es decir, se presentan con- 
figuraciones donde los grupos funcionales de los aminoácidos se alejan de la hoja de grafeno, 


siendo principalmente las cadenas laterales las que dirigen los procesos de adsorción. 
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Figura 4.1: Configuraciones de Equilibrio para cada sistema considerado, vista superior. 
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Figura 4.2: Configuraciones de Equilibrio para cada sistema considerado, vista lateral. 
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4. RESULTADOS Y DISCUSIÓN DE ADSORCIÓN DE AMINOÁCIDOS 


En los aminoácidos aromáticos como Phe y Tyr, los anillos se orientan paralelamente a 
la hoja de grafeno luego de la fisisorción, con una disposición geométrica que corresponde al 
apilamiento tipo AB Bernal [[81], bien conocido en la literatura para grafeno bicapa, donde 
los anillos se centran en un átomo de carbono de la hoja de grafeno que se encuentra debajo, 
tal y como se muestra en la Fig. (4. 11d) y (h)). Los anillos de His y de Pro no adoptan una 
configuración totalmente paralela a la hoja del grafeno (se observa una inclinación para el 
anillo de His de 17” aproximadamente, ver Fig. (4.3), esto puede deberse a la naturaleza 
de los anillos, ya que, el anillo imidazol de la histidina es un pentágono nitrogenado y el 
anillo de la prolina es una pirrolidina (estructura cíclica que está compuesta por un átomo 
de nitrógeno y cuatro de carbono, no hay presencia de dobles enlaces). En otras palabras, 
los anillos de Phe y Tyr tienen estructura hexagonal y se componen de carbono e hidrógeno, 
mientras que, His y Pro son pentágonos que incluyen nitrógenos. 

La configuración de equilibrio para Tyr e His calculada predictivamente en este trabajo, 
es corroborada posteriormente en trabajos experimentales de Krebs y colaboradores [82], 
quienes reportan una disposición paralela para el anillo de Tyr con respecto al substrato, 
mientras que para el anillo de His se reporta una inclinación. Sus estudios se llevaron a 
cabo mediante un análisis de microscopía de efecto tunel (STM), de un susbtrato de grafito 


altamente orientado (HOPG) con moléculas de His y Tyr adsorbidas. 





Figura 4.3: Configuración de Equilibrio para la molécula de Histidina. El anillo imidazol se 


encuentra inclinado 17.38". 


Los aminoácidos no aromáticos como Lys y Arg extienden sus cadenas laterales sobre 
la hoja de grafeno maximizando las interacciones tipo van der Waals, manteniendo a su vez, 
a los grupos carboxilos más apartados de la hoja de grafeno. Esto ocurre también con la 
configuración de equilibrio de la Pro, donde, el grupo carboxilo de aleja de la superficie 
de grafeno. Aminoácidos como Gly, Ala, Asp y Glu extienden sobre el grafeno sus grupos 


carboxilo y amino. 
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Figura 4.4: Energía y distancia de Adsorción para aminoácidos adsorbidos sobre grafeno. El 
driagrama de barras se divide en dos partes: superior, corresponde a los resultados de d““(Á) 
para cada aminoácido; inferior, E“ds (ev) de cada aminoácido, valores más negativos de peds. 
corresponden a sistemas más estables. Se comparan los resultados de este trabajo (barras sólidas), 
con los datos reportados por otros autores (barras con patrones). Autores de las Ref.[6] y la 
Ref.[[7] realizan estudios con hojas de grafeno dopadas con hidrógeno y considerando sólo las 
cadenas latelares (CL), respectivamente. Las distancias reportadas en la Ref. [6] corresponden a 
la distancia más cercana entre el aminoácido y el grafeno, mientras que, las distancias reportadas 
por la Ref. [7] son las distancias interplanares promedio. Las distancias reportadas en este trabajo 


corresponden a la distancia entre el centro geométrico de la molécula y el grafeno 


En el cuadro. 4.1] y, en el diagrama de barras de la Fig. se reportan los resultados 
para la energía y distancia de adsorción de los aminoácidos contemplados en este estudio. 
De acuerdo a los cálculos se obtiene la siguiente relación para las energías de adsorción: 
Arg/Grap > Lys/Grap > Tyr/Grap > His/Grap > Phe/Grap > Glu/Grap > Asp/Grap > 
Pro/Gra > Ala/Gra > Gly/Gra. Los aminoácidos más débilmente absorbidos corresponden a 
Ala y Gly mientras que los sistemas más estables corresponden a Arg/Grap y Lys/Grap. Este 
resultado está en correspondencia con los reportes experimentales de la interacción de óxido 
de grafeno con aminoácidos, Pandit et al. informó que las constantes de ligadura y los 


cambios de entalpía son mayores para los complejos de Arg y Lys (de doce aminoácidos 
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estudiados), considerándolos como los sistemas más estables [84]. 

En el cuadro. [4.1] y, en el diagrama de barras de la Fig. (4.4), se incluyen los valores 
para las energías de adsorción teóricas calculadas por otros autores. Las diferencias notables 
se deben principalmente a los métodos usados y las estructuras modeladas por cada autor. 
Por un lado, los valores reportados por Singla y colaboradores [6] basan sus estudios en la 
corrección DFT-D2 y , por el otro, Tavassoli et al. [17] estudió la interacción de las cadenas 
laterales de los aminoácidos His y Phe sobre grafeno, sin incluir en sus cálculos los grupos 
carboxilo y amina —es decir, en sus cálculos representó a los aminoácidos His y Phe me- 
diante un anillo imidazol y un tolueno, respectivamente—. Es necesario tener en cuenta, que: 
1) el método DFT-D2 considera las interacciones de todos los pares de átomos, mientras que 
DFT-D3 también considera los tríos de átomos, es decir, tienen en cuenta los efectos de tres 
cuerpos. DFI-D3 es mucho más apropiado para el estudio de interacción de moléculas con 
superficies, más aún cuando las fuerzas de largo alcance juegan un papel importante en las 
configuraciones de equilibrio y adsorción; 11) incluir los grupos carboxilo y amino modifica 
las energías de adsorción, de acuerdo a los datos del cuadro. se evidencia una subes- 
timación del 46 % y 57 % en las energías de adsorción para His y Phe. En la literatura se 
encuentran trabajos sobre los efectos en las energías de adsorción cuando se incluyen grupos 
funcionales [85]. Rochefort et al. reportó variaciones en la magnitud de las interac- 
ciones r-T y en la energía de adsorción cuando hay presencia de grupos funcionales -COOH 

De acuerdo a lo anterior, los datos reportados en este trabajo tienen en cuenta los efectos 
de tres cuerpos y además incluyen todos los átomos de los aminoácidos. Es importante men- 
cionar que en una primera etapa del presente trabajo se estudió la influencia de los grupos 
carboxilo y amino en la adsorción de His sobre grafeno, por lo que se consideró relevante 
contemplar las interacciones con los grupos funcionales de los aminoácidos [69]—se re- 
portó modificaciones en transferencia de carga, energía de adsorción, distancia de adsorción 
e incluso, modificaciones en las propiedades de transporte—. 

La distancia de adsorción entre aminoácidos y grafeno no supera los 4 Á. Las distan- 
cias reportadas en este trabajo corresponden a la distancia entre el centro geométrico de la 

Las distancias reportadas en la Ref. [6] corresponden a las distancias más cercanas entre el aminoácido y 
el grafeno, mientras que, las distancias reportadas por la Ref. [7 son las distancias interplanares promedio. Las 


distancias reportadas en este trabajo corresponden a la distancia entre el centro geométrico de la molécula y el 


grafeno. 
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4.1 Resultados de Configuración más estable 


Cuadro 4.1: Energía y distancia de Adsorción para aminoácidos adsorbidos sobre grafeno. Au- 
tores de las Ref. [6] y la Ref. [[7] realizan estudios con hojas de grafeno dopadas con hidrógeno y 


considerando sólo las cadenas latelares (CL), respectivamente 












Aminoácido poes qe Método AQ | Grafeno | Grafeno 
(eV) (A) (e) D/A S/M 
-1.49 DFI-D3 este trabajo | 0.181 S 
H; -0.63 (CL) B3LYP-D3 — — 
1S 
-0.63 DFT-D2 [6] — — 
-0.91 3.15 | DFT-D3 este trabajo | 0.106 S 
Ala -0.35 (CL) | 2.63* B3LYP-D3 — — 
Asp . 
este trabajo 
Tyr : 
este trabajo 
-2.89 3.1 | DFI-D3 este trabajo | -0.362 A M 
Arg -1.66 (CL) | 3.05* B3LYP-D3 — — — 
DFI-D3 este trabajo | 0.128 EN S 
Glu 
3.11 | DFTI-D3 este trabajo | 0.112 S 
Gly LIO" DFT-D2 [6] — — 
3.30 | DFI-D3 este trabajo | 0.100 S 
Phe 3.36* DFT-D2 [6] — — 
Pro l 
este trabajo 
Lys E 
este trabajo 


4. RESULTADOS Y DISCUSIÓN DE ADSORCIÓN DE AMINOÁCIDOS 


molécula y el grafeno —debido a que la hoja de grafeno sufre deformaciones en el proceso 
de fisisorción, para cada geometría de equilibrio, se determinó la distancia promedio en la 
dirección z de la hoja de grafeno—. Los valores calculados para distancias y energías de 
adsorción dan cuenta de la fisisorción fuerte que se lleva a cabo en los procesos de interac- 
ción entre aminoácidos y grafeno, no se observan formaciones de enlaces químicos entre los 


aminoácidos y los átomos de carbono de la hoja de grafeno. 


4.2 Densidad de Carga 


Con el objetivo de determinar si el grafeno es aceptor o donor de electrones en los procesos 
de fisisorción con aminoácidos, se realizó el análisis de carga por “Mulliken population [| En 
este trabajo se reporta una transferencia de carga entre 0.18 e y 0.1 e para el caso de grafeno 
como donor, mientras que, la transferencia de carga desde el aminoácido al grafeno tiene 
una magnitud de 0.4 e. En el cuadro. se presentan los valores de transferencia de carga 
para cada caso estudiado. Los valores reportados teóricamente, tienen el mismo orden de 
magnitud de valores medidos experimentalmente para la adsorción de moléculas 7"70/8'80 - 
tetracyano-p-quinonedimethane (INCQ) sobre grafeno (se reporta una transferencia de carga 
de 0.99 e desde el grafeno a la molécula) y, en adsorbatos más pequeños como NO» y H20 
(donde el grafeno dona 0.09 e y 0.025 e, respectivamente). La densidad de carga luego de la 


fisisorción de aminoácidos se calculó mediante la siguiente expresión: 


D= Daa/cra — Daa — Dora (4.1) 


siendo Dar /qra la densidad de carga del aminoácido adsorbido sobre grafeno; Daa la 
densidad de carga del aminoácido y, Dara la densidad de carga del grafeno. La densidad de 
carga D para cada sistema, se presenta en la Fig. (4.5). Las regiones color cyan corresponden 
a D> 0, es decir, alta densidad electrónica luego de la fisisorción; las regiones color magenta 


corresponden a D< 0, regiones de baja densidad electrónica. 


'De acuerdo a la literatura, el esquema de Mulliken constituye una opción rápida y barata (computacional- 
mente) para calcular cargas —puede ser usado para interpretar y predecir el comportamiento reactivo de una 


amplia variedad de sistemas químicos en reacciones tanto electrofílicas como nucleofílicas— 871. 


S6 


LS 





Figura 4.5: Densidad de Carga. Las regiones color cyan corresponden a D> 0, es decir, alta densidad electrónica luego de la 


fisisorción; las regiones color magenta corresponden a D< 0, regiones de baja densidad electrónica. 


e3.18) 3P PPPIsua 7'p 
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El grafeno es un donor de electrones cuando adsorbe His, Ala, Asp, Tyr, Glu, Gly, Phe 
y Pro, además de acuerdo a la redistribución de carga de la Fig. (4.5), en estos casos se 
observan regiones de alta densidad electrónica a través de los aminoácidos, mientras que, 
sobre la hoja de grafeno, se visualiza una región de pérdida electrónica. Por otro lado, el 
grafeno actúa como aceptor de electrones cuando adsorbe Arg y Lys, presentando regiones 
de alta densidad electrónica sobre la hoja de grafeno. 

En los sistemas His/Gra, Ala/Gra, Tyr/Gra, Arg/Gra, Glu/Gra, Phe/Gra, Pro/Gra y Lys/- 
Gra se forman dipolos verticales entre los adsorbatos y la hoja de grafeno —con regiones 
de baja densidad electrónica sobre la hoja de grafeno para los sistemas semi-metálicos y, 
regiones de alta densidad electrónica sobre la hoja de grafeno para sistemas metálicos—. En 
los casos particulares de aminoácidos con anillos aromáticos en sus cadenas laterales, como 
Phe y Tyr, la formación de dipolos localizados se presenta claramente debajo de sus anillos, 
dando un claro indicio de las fuerzas tipo r-7 stacking que se generan tras la interacción con 
grafeno. 

Una polarización de carga horizontal sobre la hoja de grafeno (que forma un dipolo 
longitudinal sobre la hoja de grafeno) se presenta en la adsorción de Asp y Gly, donde los 
grupos carboxilo y amino son los responsables de los dipolos transversales que rompen la 
simetría de dipolo vertical de los demás aminoácidos absorbidos. 

Los electrones r delocalizados del grafeno son los que participan en la re-distribución 
electrónica de los sistemas. Experimentalmente, la formación de dipolos en la interacción 
entre aminoácidos y grafeno ha sido reportado por Mallineli et al. —en sus trabajos se reporta 
cargas parciales negativas (o positivas) sobre (o debajo) de una hoja de grafeno en los sitios 


de adsorción de los aminoácidos— [881. 


4.3 Densidad de Estados 


Con el fin de obtener información de las propiedades físicas de los sistemas estudiados, se 
calcularon las densidad de estados totales para el grafeno antes y después de la adsorción, 
ver Fig. y Fig. (4.7), respectivamente. 

Antes de la adsorción de aminoácidos, el grafeno presenta una densidad de estados que 
corresponde a un semimetal —sin estados en el nivel de Fermi—, donde los orbitales p, de 
los átomos de carbono, son los que aportan estados en energías entre -5 eV y 5 eV, ver la 
Fig. línea verde. Los estados entre -3 eV y -10 eV son contribuciones de orbitales p.,,, 
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Total Grafeno 
HH Orbital s 
— Orbital p_ + p, 


-— Orbital p, 


DOS (estados Vol ' eV”) 





Figura 4.6: Densidad de Estados para el grafeno antes del proceso de adsorción. La densidad de 
estados total del grafeno se representa por el color gris sólido. La contribución de los estados p., 


Py + Pu y 5, se representan por la línea verde, línea negra y azul sólido, respectivamente. 


Py Y Pz (Py + Pa línea negra en la Fig. (4.6)). A bajas energías, entre los -20 eV y -10 eV, 
predominan los estados de los orbitales s. 

Luego de la adsorción, la interacción con los aminoácidos induce la aparición de estados 
cerca a la energía de Fermi del grafeno en la banda de valencia —entre 0 eV y 2.2 eV—, 
los picos en energías de -1.2 eV, -1.8 eV, -1.3 eV, -2.11 eV y -1.49 eV provienen de His, 
Asp, Tyr, Phe y Pro, respectivamente. En el caso de Ala, Glu y Gly los estados más cercanos 
al nivel de Fermi están a -2.3 eV, los cuales se superponen al primer pico de la banda de 
valencia de la DOS del grafeno cuando no ha adsorbido moléculas. 

El carácter semimetálico del grafeno se modifica cuando adsorbe Lys y Arg. En este caso, 
el grafeno se metaliza, ya que las moléculas transfieren cargas de -0.464 e y -0.362 e, para 
Lys y Arg, respectivamente. Los demás aminoácidos considerados His, Ala, Asp, Tyr, Gly, 
Glu, Phe y Pro, mantienen el carácter semimetálico del grafeno. 

Un aspecto relevante, consistió en identificar si los estados cercanos a la E de la banda 
de valencia provenían de los grupos funcionales -COOH y -N Hb Ó, de las cadenas laterales de 


las moléculas —lo que podría dar cuenta del papel que desempeñan los grupos funcionales 
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Graphene 
aa His/Grap 
—  Ala/Grap 
— — Asp/Grap 

-+ Tyr/Grap 


«—= Phe/Grap 
E «—= Pro/Grap 
AN | Lys/Grap 


DOS (estados Vol eV”) 





Figura 4.7: Densidad de estado de los sistemas luego de llevarse a cabo el proceso de adsorción. 
La densidad de estados total del grafeno antes de la adsorción (grafeno puro), se representa por 


el color gris sólido. 


en los procesos de adsorción—. Para ello se calcularon las densidades de estados parciales 
de los grupos funcionales y las cadena laterales para cada aminoácido fisisorbido, ver la 
Fig. (4.8). 

De acuerdo a los cálculos, los picos cercanos a la E provienen principalmente de la 
interacción con los grupos funcionales, ya que sólo dos de los diez aminoácidos estudiados 
presentan picos atribuidos únicamente a sus cadenas laterales. Según la Fig. los estados 
provienen de 1) cadenas laterales, para His y Tyr, 11) -COOH y -NH») para Asp y Ala, 111) 
-NH) para Arg, Lys, y Gly, 1v) cadena lateral y - NH» Glu, Phe y Pro. 

Los grupos funcionales juegan un papel preponderante en la aparición los estados cerca 
al nivel de Fermi, no sólo porque en la mayoría de las moléculas estudiadas hay participación 
directa —ocho de los diez aminoácidos estudiados introducen estados cerca al nivel de fermi 
provenientes de los grupos funcionales—, sino porque, tal y como se reportó en los estudios 


titulados “Computational study of transport properties of graphene upon adsorption of an 
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amino acid: importance of including -NHy and -COOH groups”[69]] y “Theoretical study of 
the adsorption of histidine amino acid on graphene”[68], la inclusión de los grupos funcio- 
nales 1) aumenta la energía de adsorción y la transferencia de carga y, 11) desplaza hacia la 
energía de Fermi los estados que provienen de la cadena lateral. 

En los artículos referenciados en el párrafo anterior, se sintetizan los resultados de la 
adsorción de histidina con y sin grupos funcionales, ya que para esa fecha, los primeros tra- 
bajos de adsorción de aminoácidos sobre grafeno sólo consideraban los efectos de la cadena 
lateral (que en este caso específico para histidina, corresponde a un anillo imidazol). Nues- 
tro modelo muestra diferencias significativas para energía de enlace y transferencia de carga 
entre el aminoácido con y sin grupos carboxilo y amino, además del desplazamiento en 0.29 


eV hacia Ep del pico de la cadena lateral una vez se incluyen los grupos funcionales. 
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PDOS (estados eV *? Vol 1) 
PDOS (estados eV ** Vol -?) 





E (eV) 


Figura 4.8: Densidad parcial de estados (PDOS). La densidad de estados total para cada sistema 
AA/Grap, se representa con color cyan. Se reporta la contribución a la DOS de las cadenas 
laterales (CL con línea verde), los grupos carboxilo (línea color magenta) y grupos amino (color 
azul sólido). Se incluye en cada gráfica la DOS del grafeno antes de la adsorción (línea sólida 


negra) 
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Efecto de la adsorción de 
aminoácidos en el transporte 
electrónico de una hoja de 
srafeno 


Con el fin de determinar los efectos de la adsorción de aminoácidos en las propiedades 
de transporte electrónico de una hoja de grafeno, se calcularon los espectros de transmisión y 
las curvas de corriente-voltaje para once dispositivos: Grap, His/Grap, Ala/Grap, Asp/Grap, 
Tyr/Grap, Arg/Grap, Glu/Grap, Gly/Grap, Phe/Grap, Pro/Grap y Lys/Grap (ver, Fig. (5.1). 
Tal y como se mencionó en el capítulo 3, las regiones centrales (C) de los dispositivos mo- 
delados de la Fig. (5.1), corresponden a las configuraciones de equilibrio presentadas en el 
capítulo 4. A continuación discutimos los resultados principales del estudio de transporte. 

La corriente que fluye a través de cada sistema modelado, se determina mediante la inte- 
gración de la transmisión electrónica sobre una ventana de tensión. En este orden de ideas, el 
aumento, disminusión o supresión total de picos en las curvas de transmisión, brindará infor- 
mación sobre el incremento o el decrecimiento de la corriente. En la Fig. se presentan 
los espectros de transmisión para los voltajes de polarización de -0.8 V, -1.5 V, -2 V, 0.8 V, 
1.5 V y 2 V. 
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(f) 


(b) 


(h) 


(c) 


(1) 





0) 





Figura 5.1: Dispositivos modelados (vista superior). La región C de cada dispositivo se deter- 


minó en el capítullo 4. 


De acuerdo a los cálculos, la transmisión depende del voltaje de polarización (6 V;) 
que se implementa para modelar cada dispositivo, ya que, una variación en la magnitud o 
en la polarización (positiva/negativa) del V;,, conlleva a que algunos picos de la curva de 
transmisión de electrones se atenúen, se supriman o aumenten. En el sistema Phe/Grap, por 


ejemplo, la integración de la transmisión para V; positivos es menor que para V; negativos. 
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Figura 5.2: Transmisión de los sitemas modelados para voltajes bias de -0.8 V, -1.5 V, -2 V, 0.8 V, 1.5 V y 2 V. 


5. EFECTO DE LA ADSORCIÓN DE AMINOÁCIDOS EN EL TRANSPORTE 
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Analizando específicamente las curvas para voltajes de polarización negativos para Phe 
(ver, Fig. (5.2), la integración de la transmisión a V, = —0,8 V es menor que para V, = —2 
V, donde el pico A (ver, Fig. para Phe), aumenta cuando el voltaje aumenta, lo que 
implica que, bajo un voltaje de polarización de -2 V los electrones con esa energía particular 
van a tener mayor probabilidad de transmitirse del electrodo izquierdo al derecho. 

A diferencia de Phe/Grap, en el sistema Lys/Grap la integración de la transmisión para 
V; positivos es mayor que para V;, negativos. La dependencia de la transmisión respecto 
al V, se observó también en trabajos sobre sensado de ADN, desarrollados por Souza y 
colaboradores [89] [90]. 

En las Fig. y Fig. (5.441) se presentan las curvas de corriente vs voltaje bias (1-V;) 
para el grafeno antes y después de la adsorción de cada aminoácido. En la primera figura, 
se muestran los valores de corriente para voltajes bias entre O V y 1 V, mientras que, en la 


segunda se graficaron los datos para voltajes bias entre -2V y 2V. 


p-a His/Grap 
Ala/Grap 
Asp/Grap 
Tyr/Grap 


Arg/Grap 
Glu/Grap 
Gly/Grap 
Phe/Grap 


Pro/Grap 
Lys/Grap 
Grafeno 





Figura 5.3: Curva de corriente-voltaje para los dispositivos modelados (V; entre O V y 1 V). 


La curva de corriente para el grafeno antes de la adsorción (línea negra sólida en las 
Fig. (5.3) y Fig. (5.44)) corresponde a una curva típica de un dispositivo semiconductor de 
gap cero, exhibiendo un comportamiento no lineal y simétrico para voltaje bias positivos y 


negativos —simétrico con respecto a la magnitud de la corriente bajo inversión de voltaje 
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bias—. Una relación no lineal entre la corriente y el voltaje bias aplicado a una hoja de 


erafeno, ha sido reportado también por varios autores 1791. 


a-—a His/Grap 
. —= Ala/Grap 
A--A Asp/Grap 
Rx Tyr/Grap 
.e- + Arg/Grap 
—  Glu/Grap 
— Gly/Grap 
o- O Phe/Grap 
—— Pro/Grap 


e—e lys/Grap 
e—o Grafeno 





0 
V,(V) 


Figura 5.4: (a) Curva de corriente-voltaje para los dispositivos modelados (V; entre -2 V y 2 V). 
(b) diferencia de corriente entre el grafeno antes y después de la adsorción con cada aminoácido 


Al = IGrap == TAA/Gra: 


Una vez que se lleva a cabo la adsorción de cada amino ácido, el grafeno presenta cam- 
bios significativos en las curvas I-V, producto de la interacción con las moléculas sobre su 
superficie, donde los orbitales moleculares pueden permitir o inhibir el transporte electrón1- 
co a través de la hoja de grafeno. De acuerdo a la Fig. (5.3), los valores de corriente son 
pequeños cuando el voltaje bias se encuentra entre -0.6 V y 0.6 V —para nueve de los diez 


aminoácidos estudiados la corriente calculada no supera los 0.04 yA—. Para este mismo 
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rango de voltaje de polarización, el sistema Lys/Grap registra un valor de corriente levemen- 
te mayor con respecto a los demás sistemas, de 0.1 yA. En la Fig. (5.4p) se observa que para 
magnitudes de V;, superiores a los 0.6 V la corriente calculada depende de cada aminoáci- 
do. Adicionalmente, hay una asimetría en las curvas de corriente respecto a la inversión 
del voltaje de polarización, esto se debe a las configuraciones de equilibrio que adoptan los 
aminoácidos al interactuar con el grafeno y la asimetría intrínseca de su densidad de carga. 

La asimetría de las curvas I-V son de gran interés en el campo de la electrónica molecular 
y el estudio de nanodispositivos. Una molécula asimétrica tendrá orbitales moleculares con 
una distribución de carga asimétrica. Esto en la práctica induce un acoplamiento diferente 
con los electrodos que puede conducir a un perfil de voltaje asimétrico [91]. Así mismo, 
hay moléculas que presentan características propias de un diodo rectificador, es decir, que 
bloquean el transporte de carga en una dirección y permiten el flujo en la dirección opuesta, 
un ejemplo son las monocapas de moléculas aromáticas ensambladas (SAM) en grafeno que 
mostraron un comportamiento de diodo molecular 1931. 

Con el fin de especificar las variaciones de corriente de cada sistema respecto al grafeno 
sin molécula, se graficó la diferencia de corriente entre el grafeno antes y después de la 


adsorción con cada aminoácido, se definió la diferencia de corriente (Af) como: 


Al = larap e TAA/Gra (5.1) 


Donde lfarap es la corriente del grafeno antes de la adsorción y L1A/Gra es la corriente 
del grafeno después de adsorber un aminoácido. 

En la Fig. (5.4p) se observa que los aminoácidos absorbidos sobre grafeno como la Lisina 
(Lys) pueden actuar como un diodo rectificador, permitiendo mayor flujo de carga para V, 
positivos y bloqueando para V; negativos, se calculó una diferencia de corriente de 2.19 
HA y 0.12 y A para voltajes de polarización de 2 V y -2 V, respectivamente. Moléculas 
como fenilalanina (Phe) y ácido aspártico (Asp) presentan un comportamiento inverso al 
caso de Lys, ya que, permiten mayor paso de corriente para valores negativos de voltaje de 
polarización e inhiben el paso de corriente para valores positivos —la magnitud de A/ es 
de 0.81 yA y 2.74 A para voltajes de polarización de 2 V y -2 V, respectivamente—. Los 
sistemas con perfiles de polarización más simétricos respecto al voltaje de polarización en 


comparación al resto de aminoácidos son Pro/Grap y Arg/Grap. 
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De acuerdo a los datos obtenidos, el orden de magnitud de la diferencia de corriente 
para todos los casos no es menor a los yA, siendo consistente con datos experimentales de 
dispositivos tipo FET basados en grafeno. En el cuadro. podemos referenciar algunas 


diferencias de potencial medidos para la detección de bacterias y otras biomoléculas. 


Cuadro 5.1: Superior: Valores de la diferencia de corriente A/ de los dispositivos modelados en 


este trabajo. Inferior: Valores de la diferencia de corriente registrados en dispositivos tipo GFET 


experimentales. 
Dispositivo modelado en este trabajo 
A KA) A KA) A KA) Aminoácido 
(V¿=-2 VW) (V,=-1.5 V) (V¿=-1 V) 
0.32 0.83 0.28 His 
0.83 1.03 0.34 Ala 
1.94 1.54 0.69 Asp 
-0,47 0.27 0.17 Tyr 
-0.56 0.23 0.08 Arg 
-0,.53 0.47 0.17 Gly 
0.39 0.44 0.13 Glu 
2.67 2.15 0.93 Phe 
-0.66 0.05 0.03 Pro 
0.11 0.31 0.15 Lys 


Trabajos Experimentales. Dispositivos basados en grafeno 


Referencias  ATKuA) Biomolécula  TipodeMolécula 
5 E.Coli Bacteria 
4 ADN ácido nucleico 
[22] 10 ADN 
[97] 1 ADN 


Un aspecto interesante de analizar, hace referencia a la sensibilidad del dispositivo mo- 
delado. En este trabajo, el mecanismo de detección de aminoácido se basa en cambios de 
corriente, por lo que, identificar los porcentajes de incremento/decrecimiento de la señal de 


corriente de grafeno una vez adsorbe cada aminoácido, resulta apropiado para dar cuen- 
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ta de sensibilidad del dispositivo —tal y como lo proponen varios autores, en estudios de 
transporte de grafeno aplicados al sensado de ADN [9Q]—. Se define el porcenaje de 


incremento/decrecimiento de la señal de grafeno, como: 


Al 


La 


API= 


100 (5.2) 





donde Y% PI es el porcentaje de incremento de la señal de grafeno, A/ se definió en la 
ec. y [Gray es la corriente del grafeno antes de la adsorción. En la Fig. se presentan 
los porcentajes de incremento de la señal para voltajes de polarización de 1 V, 1.5 V y 2 V. 
El % PI depende del voltaje bias que se esté considerando, así por ejemplo, a 1 V, el sistema 
Lys/Grap, presenta un incremento del 257 % de la señal de corriente respecto al grafeno puro, 
mientras que, el incremento a 1.5 V y 2 V es de 133 % y 90 %, respectivamente —en otras 
palabras, el dispositivo para Lys/Grap es más sensible a 1 V que a 1.5 V y 2 V— , 

Para un voltaje de polarización de 1 V, los sistemas que presentan un incremento de la 
señal mayor al 75 % son Lys/Grap, Asp/Grap y Ala/Grap con porcentajes de 257 %, 110% 
y 85 %, respectivamente. Por otro lado, para un voltaje extremo de 2 V, el porcentaje de 
incremento para Glu/Gra y Gly/Grap es de aproximadamente del 3 %, lo que convierte a la 
Glu y Gly en los aminoácido más difícil de distinguir (en comparación al grafeno solo) para 
ese voltaje. 

Las curvas I-V calculadas para el dispositivo propuesto, dan cuenta de las potencialidades 
del grafeno en la construcción de mecanismos de sensado selectivos y sensibles a pequeñas 
variaciones de carga locales, ya que, por un lado, se observan diferentes curvas de corriente 
para cada aminoácido, y por el otro, el orden de magnitud de A/ corresponde a dispositivos 


experimentales. 
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Figura 5.5: Porcentaje de incremento/decrecimiento de la señal de grafeno Y PI, para voltajes 
de polarización de -1.5 V, -1 V y 1 V. 
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Conclusiones 


En este trabajo se estudió teóricamente mediante cálculos DFT y NEGE, la interacción de 
diez aminoácidos con una hoja de grafeno. El trabajo se dividió en dos partes 1) la adsorción 
de aminoácidos sobre grafeno (His, Ala, Asp, Tyr, Arg, Gly, Glu, Phe, Pro y Lys) y, 11) estudio 
del transporte electrónico del grafeno antes y después de la adsorción con cada aminoácido, 
donde se calcularon curvas de transmisión y corriente-voltaje. 


Sobre el estudio de adsorción 


e De acuerdo a los resultados, los diez aminoácidos se fisisorben sobre la superficie de 
erafeno, es decir, en los procesos de interacción entre substrato y adsorbato, no se 


generan enlaces covalentes. 


e Las distancias de adsorción no superan los 4 Á—-las distancia entre substrato y adsor- 
bato para cada sistema, se determinó como la distancia entre el centro geométrico de 


la molécula en z y la hoja de grafeno—. 


e Las energías de adsorción presentan el siguiente orden de jerarquía: Er 1arap > ELys/Grap > 


ads ads ads ads ads ads ads 
E ryr/Grap > Crap > e rhlGras > Ectu/Grap > can ES rro iCra > Ea lCrá A 
Edy /ara: LOS Cálculos de energía de adsorción concuerdan con resultados experimen- 
tales de interacción de aminoácidos y grafeno, donde los sistemas más estables corres- 


ponden a Arg/Grap y Lys/Grap. 
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6. CONCLUSIONES 


e En cuanto a las configuraciones de equilibrio, los aminoácidos se fisisorben sobre gra- 
feno sin un lugar privilegiado para los grupos funcionales, es decir, hay moléculas 
donde los grupos funcionales se enfrentan a la hoja de grafeno (como en el caso de 
Tyr, Gly, His ), mientras que en otras moléculas, los grupos funcionales se alejan (co- 
mo en el caso de la Pro, Lys y Phe), regidos básicamente por la transferencia de carga 


y los dipolos asociados. 


e Los aminoácidos con cadenas laterales aromáticas como His, Tyr y Phe, localizan sus 
anillos paralelamente a la hoja de grafeno (excepto para el anillo imidazol de His, que 


presenta un ángulo de inclinación 17.38". 


e En todos los procesos de fisisorción se calcula una transferencia de carga desde o hacia 
la hoja de grafeno, con la conformación de dipolos locales longitudinales y transver- 
sales entre la moléculas y los enlaces p del grafeno. Experimentalmente, resultados 
de interacción de aminoácidos con grafeno muestran evidencia de la conformación de 


dipolos de carga. 


e El grafeno es un donor débil de electrones cuando adsorbe His, Ala, Asp, Tyr, Gly, 
Glu, Phe y Pro, con transferencias de carga hacia las moléculas que no superan los 0.2 
e. En el caso específico de Lys y Arg se presenta una transferencia de carga desde las 
moléculas hacia el grafeno de -0.464 e y -0.362 e, respectivamente. La carga ganada 
por la hoja de grafeno, se redistribuye sobre la superficie del material, contribuyendo 


a la formación de dipolos de carga. 


e De acuerdo a las densidades de estado totales y locales, los procesos de fisisorción 
introducen estados cerca al nivel de fermi del grafeno. Los estados provienen de los 
grupos funcionales de los aminoácidos y no sólo de las cadenas laterales de los mis- 


mos, por lo que toda modelización de éstas moléculas debe incluirlos. 


e Los estados de la banda de valencia más cercanos al nivel de Fermi se encuentran en 
energías de -1.2 eV, -1.8 eV, -1.3 eV, -2.11 eV y -1.49 eV para His, Asp, Tyr, Phe y 
Pro, respectivamente. En el caso de Ala, Glu y Gly los estados más cercanos al nivel 
de Fermi están a -2.3 eV, los cuales se superponen al primer pico de la densidad de 


estados del grafeno cuando no ha adsorbido moléculas. 
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En todos los casos, la contribución en la densidad de estados parciales de los gru- 
pos carboxilos, se presenta mediante la aparición de dos picos separados 1.2 eV en 
la banda de valencia (los dos primeros picos más cercanos a la banda de valencia). 


Específicamente, el pico en -3.8 eV en la PDOS proviene de los grupos carboxilos. 
La adsorción de Arg y Lys presenta picos en -3.1 eV provenientes de los grupos amino. 


Los grupos funcionales juegan un papel preponderante en la aparición de los estados 
cerca al nivel de Fermi, no sólo porque en la mayoría de las moléculas estudiadas hay 
participación directa —ocho de los diez aminoácidos estudiados introducen estados 
cerca al nivel de fermi provenientes de los grupos funcionales—, sino porque, tal y 
como se reportó en un estudio titulado Computational study of transport properties of 
graphene upon adsorption of an amino acid: importance of including "NH and —-CO- 
OH groups publicado en el año 2016, la inclusión de los grupos funcionales aumenta 


la energía de adsorción y la transferencia de carga. 


Las densidades de estado muestran la metalización del grafeno luego de la adsorción 
de Lys y Arg, producto de la transferencia de carga de las moléculas al grafeno. Para 
el resto de los aminoácidos considerados His, Ala, Asp, Tyr, Gly, Glu, Phe y Pro, se 


mantiene el carácter semimetálico del grafeno. 


Sobre el estudio de transporte 


Se calcularon las curvas de transmisión y corriente-voltaje del dispositivo tipo GFET 
con electrodos de grafeno. Los perfiles de las curvas I-V son diferentes para cada 
aminoácido, lo que deja en evidencia la potencialidad del material en la construcción 


de dispositivos selectivos de sensado. 


Las curvas de transmisión se modifican para cada aminoácido dependiendo del poten- 
cial de polarización. Aminoácidos como Phe y Asp tienen una mayor transmisión para 
voltajes de polarización negativa, mientras que la Lys permite mayor transmisión para 


voltajes positivos. 


Las diferencias en las curvas de transmisión electrónica para valores de polarización 
positivos y negativos, son producto de la densidad de carga asimétrica propia de la 
molécula (la cual se acopla de manera diferente con cada electrodo, mejorando o dis- 


minuyendo la transmisión electrónica de un electrodo a otro). 
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6. CONCLUSIONES 


e Se calcularon las diferencias de corriente, entre el grafeno puro y el sistema aminoáci- 
do adsorbido por grafeno. Las diferencias obtenidas en este trabajo, son del mismo 
orden de magnitud de las medidas experimentales para dispositivos GFET —se han 
reportado dispositivos sensibles a variaciones de corriente de alrededor de 0.05 yA—., 
Para un voltaje de polarización de 1.5 V, la menor diferencia de corriente para el dis- 
positivo propuesto es de 0.22 yA y el máximo de 1.64 A, lo cual lo convierte en 


perfectamente funcional. 


e La asimetría en los perfiles de las curvas de corriente-voltaje develan otras posibles 
aplicaciones de aminoácidos en la construcción de nanodispositivos electrónicos ba- 
sados en grafeno, como por ejemplo, desarrollo de diodos rectificadores moleculares. 
Los aminoácidos Phe y Lys, presentan un comportamiento claro de diodo rectificador, 
en el primer caso, se inhibe el paso de corriente para voltajes de polarización positi- 
vos, mientras que se permite en valores negativos. En el segundo caso, se presenta el 


proceso Inverso. 


Perspectivas 
Para trabajo futuros, resultaría interesante explorar las propiedades de transporte de siste- 
mas GFET análogos al propuesto en esta tesis. Materiales como el óxido de grafeno pueden 


llegar a presentar otros comportamientos cuando adsorben aminoácidos. 
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Autoensamblaje molecular sobre 
srafeno 


En el marco del desarrollo del último año de doctorado, se inició un trabajo de adsor- 
ción de moléculas de Di-carbonitrilo-polifenilo sobre grafeno. Específicamente, el estudio 
de moléculas autoensambladas. Este trabajo se realizó en colaboración con el investigador 
Juan Carlos Moreno, quién suministró información de los datos experimentales||| 

Existen diferentes caminos para controlar dispositivos basados en grafeno: determinando 
cambios en corriente y conductancia, manipulando el spin ó, estudiando la modificación de 
los valles del espectro de energía de la banda de conducción (esta ultimo llamado valletrónica 
ó “valeytronic” en inglés) [99]. En el proceso de controlar un dispositivo, dependiendo 
de la aplicación requerida, la funcionalización del grafeno juega un papel preponderante. En 
este sentido el autoensamblaje molecular es una estrategia prometedora, ya que, existe una 
gran variedad de moléculas que se pueden adsorber sobre grafeno —i1ncluso átomos aislados 
con alto momento magnético—, dando lugar a la exploración de materiales con aplicaciones 


en la espintrónica. 


Es necerario aclarar, que presentamos los resultados parciales de los cálculos teóricos, no es un trabajo 
finalizado. 
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A. AUTOENSAMBLAJE MOLECULAR SOBRE GRAFENO 


A.1 Introducción 


La funcionalización de grafeno con moléculas orgánicas ha generado gran interés en el desa- 
rrollo de dispositivos híbridos con aplicaciones en electrónica y optoelectrónica. En la ac- 
tualidad, una estrategia para lograr dicha funcionalización, consiste en formar monocapas 
autoensambladas de moléculas orgánicas (SAM por sus siglas en inglés Self-assembled mo- 
nolayers). 

El autoensamblaje molecular es el proceso de asociación espontánea de moléculas (ó 
partes de moléculas) para la formación de estructuras más complejas como supramoléculas, 
estructuras porosas, entre otras. La variedad de moléculas orgánicas y metales utilizados en 
los procesos de autoensamblado sobre superficies, permite diseñar a “priori” nanoestructuras 
de poros con dimensiones y geometrías específicas. 

De acuerdo a la literatura, entre los conceptos claves para estudiar SAMs sobre diferentes 


substratos se tienen: 


i. Componentes. Hace referencia al grupo de moléculas o segmentos (también conocidos 
como ligandos) que interactúan entre ellas para conformar la SAMs. Los ligandos 
pueden ser iguales o diferentes, los cuales, por lo general tienen grupos funcionales con 
mayor afinidad al substrato. Entre los grupos funcionales típicos se tienen los hidroxilo 
(-OH), amino (-NH>), carboxilo (-COOH) o carbonitrilo (-CN). Cabe mencionar aquí, 


que hay SAMs que requieren átomos metálicos individuales. 


ll. Interacciones. El autoensamblado ocurre cuando las moléculas interactúan entre sí a 
través de un equilibrio de fuerzas atractivas y repulsivas. Estas interacciones son ge- 
neralmente débiles (es decir, comparables con las energías térmicas) y no covalentes 
(interacciones de van der Waals y Coulomb, interacciones hidrofóbicas y enlaces de 
hidrógeno), aunque también, se pueden contemplar enlaces covalentes relativamente 
débiles (enlaces de coordinación). Las interacciones se dan entre los ligandos y/o entre 
ligando y substrato. En el caso de SAMs con átomos metálicos individuales, se gene- 


rarán otras interacciones, tales como, metal-ligando, metal-metal y metal-substrato. 


En el último año de doctorado, se llevó a cabo el modelado computacional mediante 
la teoría del funcional de la densidad (DET), de la formación de redes autoensambladas con 
ligandos de Di-carbonitrilo-polifenilo NC-Ph¿-CN (P6), NC-Ph3-CN (P3) sobre un substrato 
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A.2 Proceso experimental de elaboración de las SAMs 


de grafeno (ver ligandos en la Fig. (A.1)). Se estudiaron sistemas con y sin interacción metal- 


ligando. El metal que se usó en el proceso de SAMSs fue cobalto. 


P6 
ha pa es ss ss ps 
A .. .—. 
pr .. e. .. e. ? , 





Nitrilos 


P3 
pipi — 


Figura A.1: Ligandos usados en el modelado de autoensamblaje molecular sobre grafeno. Áto- 
mos de carbono (C), hidrógeno (H) y nitrógeno (N) en color café, gris y azul, respectivamente. 


En la imágen se destacan los grupos nitrilos, importantes en el proceso de autoensamblado. 


A.2 Proceso experimental de elaboración de las SAMs 


Las SAMs modeladas, fueron sintetizadas experimentalmente por el instituto de materiales 
avanzados de la Universidad de Groningen, Holanda]|'| De acuerdo a trabajos experimentales 
y teóricos, las moléculas depositadas sobre superficies poco reactivas, como grafeno, da 
lugar a SAMSs cuyas dimensiones de poro y gran estabilidad dependen principalmente de 
la interacción entre ligandos, ya que son relativamente fuertes y altamente direccionadas 
(abarcando interacciones no covalentes, metal-ligando y enlaces de hidrógeno). 

En esta sección, se describe brevemente el proceso de síntesis experimental de las molécu- 


las autoensambladas contempladas en este estudio. 


1. Preparación del sustrato: se realiza un recocido (annealing) del sustrato a una tempe- 
ratura de 520 K (este annealing es necesario para remover agua u otras moléculas que 


podrían haberse adsorbido al grafeno durante el contacto con el aire). 


! donde el investigador Juan Carlos Moreno y Meike Stóúhr se contactaron para llevar a cabo un trabajo 


conjunto 
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A. AUTOENSAMBLAJE MOLECULAR SOBRE GRAFENO 


2. Sublimación de moléculas: los ligandos (P6 o P3) se depositaron en el sustrato de gra- 
feno sobre cobre en la dirección (111) en ultra alto vacío (UHV) a temperatura ambien- 
te. Las moléculas se depositaron por sublimación desde una celda Knudsen calibrada 
para llevar a cabo una cobertura de una monocapa. Las imágenes de microscopía de 
efecto túnel de barrido (STM, por sus siglas en inglés scanning tunneling microscopy) 


de las estructuras autoensambladas se muestran en la Fig. (A.2). 


En el caso de las SAMSs que contemplan átomos individuales, se tiene una proceso poste- 
rior adicional: las superficies con las moléculas adsorbidas se exponen por 10 segundos a un 
rayo de átomos metálicos (átomos de cobalto Co en este caso) generados por un evaporador 
(todo el proceso se lleva a temperatura ambiente). Las imágenes que STM de las estructuras 
autoensambladas se muestran en la Fig. (A.2). 


A.3 Objetivos 


Los objetivos se determinaron de acuerdo a cada sistema, ya que se estudiaron estructuras 


autoensambladas con y sin átomos de Co. 


e Sistemas sin interacción metal-ligando (sin Co). Se contemplaron SAMs con ligandos 
de P3 y P6 por separado. Los objetivos giraron en torno a 1) Determinar la distancia 
de adsorción entre grafeno y moléculas P3. 11) Determinar la orientación teórica de los 


anillos de la molécula P6 con respecto a la hoja de grafeno. 


e Sistemas con interacción metal-ligando (con Co). Se contemplaron SAMs con ligandos 
de P3 y P6 por separado. 1) Determinar la energía total de diferentes redes e identificar 


configuraciones de mínima energía. 


A.3.1 Detalles computacionales 


El estudio ab initio se enmarca en cálculos DFT, implementando pseudopotenciales con- 
servadores de la norma del software Openmx3.8. El potencial de intercambio y correlación 
usado fue GGA-PBE con DFT-D3 (que incluye interacciones van der Waals) y, el valor de 
la energía de corrección para el potencial de Coulomb (U) de los estados d del Co es de U=4 


eV. El criterio de convergencia para la relajación fue de 0.0004Har/Bohr. 


SU 


A.3 Objetivos 





Figura A.2: Redes experimentales autoensambladas. (a) P3 autoensamblado sobre grafeno; (b) 


P6 autoensamblado sobre grafeno; (c) P3 y Co sobre grafeno; P6 y Co sobre grafeno. 


Para calcular la energía y distancia de adsorción —del sistema P3 autoensalblado sobre 
grafeno— se llevó a cabo dos relajaciones: 1) Relajación del substrato, átomos de la hoja de 
grafeno sin molécula, 11) Relajación parcial, las moléculas se dispusieron a diferentes alturas 
desde la hoja de grafeno —previamente relajada— y, para cada distancia considerada, se 
efectúo una segunda relajación parcial, donde sólo los átomos de las moléculas eran libres 
de moverse. 


La energía de interacción se calculó mediante la expresión: 


pint _— pot o pref (A.1) 


Donde E** es la energía total del sistema para una separación substrato-adsorbato h. E”% es 
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la energía total cuando no hay interacción entre substrato-adsorbato. La energía de adsorción 


E“4s será el mínimo valor de la energía de interacción. 


A.4 Resultados parciales sistemas sin cobalto 


Se determinaron las configuraciones de equilibrio para las moléculas P6 y P3 sobre grafeno. 
En la Fig. (A.3|se presentan las celdas unitarias para los sistemas considerados (configura- 


ciones de mínima energía) y en el cuadro. [A.1]se especifican las dimensiones de la celdas. 





Figura A.3: Sistemas sin Co. (a) P3 sobre grafeno: estructura tipo chevron; (b) P6 sobre grafeno: 


compacta hexagonal. El substrato de grafeno en los dos casos se ha respresentado en color verde. 


Luego de la relajación, para el sistema P3 los anillos fenilos no se giran y los grupos 


nitrilos se orientan hacia los hidrógenos de moléculas vecinas formando una red tipo chevron. 
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A.4 Resultados parciales sistemas sin cobalto 


La distancia de adsorción entre las moléculas y el grafeno fue de 3.26 Á, valores acordes a 


los reportados experimentalmente. 


Para el sistema P6, el autoensamblado presentan un acoplamiento dipolar paralelo de los 
nitrilos de tres moléculas consecutivas (estructura llamada compacta hexagonal, donde la 
mayor carga electrónica para el enlace -CN se corresponde al N), en este caso, se presenta 
un cambio de acoplamiento entre grupos nitrilos que coincide con los desplazamientos de 
cada paquete de tres moléculas P6. Para P6, los anillos fenilos se giran unos con respecto a 
otros —20* con respecto al grafeno—, ver Fig. Experimentalmente, los anillos fenilos 
presentan una organización similar a la encontrada en fase gas (orientaciones diferentes para 
cada anillo fenilo), de tal forma que, hay un desacople entre el substrato y la moléculas 
(además no hay hibridización con el substrato). Teóricamente, si la distancia entre adsorbato 
y substrato disminuye —se forza el sistema—, las moléculas tienden a tener una disposición 


totalmente paralela con respecto a la hoja de grafeno. 


Ge qa 
¿7 Y DS 


dl 


Lor y 
uy WS Y 





Figura A.4: P6 sobre grafeno.Arriba: vista frontal de las moléculas, se observa claramente la 
rotación de los anillos. Abajo: esquema del acoplamiento dipolar paralelo (resaltado en rojo) y 


antiparalelo (resaltado en gris). 


En el cuadro.|A.1]se presentan las distancias de adsorción. 
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Cuadro A.1: Sistemas sin Co. 0 es el ángulo entre los vectores a y b. 


Molécula | Dimensiones Celda | Distancia de adsorción | Orientación de anillos 
(A) (A) 


P3 ld] =17,04 0 
B| = 13,68 
ld = 25,00 


O =98,9 
P6 ld] = 20,48 | 209 
b| = 28,40 
ld = 25,00 
O =82,41 





A.53 Resultados parciales sistemas con cobalto 


Se estudiaron cuatro configuraciones moleculares de ligandos y Co en fase gas, y se discutió 
en términos de energía total y densidad de carga la posibilidad de su síntesis experimental. 


En la Fig. (A.5|se presentan las cuatro configuraciones, y en el cuadro. [A.2]se presentan 
el paraámetro de las celdas cuadradas. 


Cuadro A.2: Sistemas con Cobalto en fase gas. 


Sistema Disposición de Ligandos 


P3-Cuadrado cuadrada 


P3-Inter intercalada 
P6-Cuadrado cuadrada 


P6-Inter intercalada 





P3-Cuadrado: Un átomo de Co central ligado a cuatro moléculas P3. En este complejo 
no hay interacción entre ligandos. El autoensamblaje da lugar a la formación de una red 


cuadrada plana con átomos de Co en los vértices. 
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A.5 Resultados parciales sistemas con cobalto 
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Figura A.5: Sistemas con Co (a) P3 sobre grafeno disposición cuadrada; (b) P3 sobre grafeno 
disposición entrelazada; (c) P6 sobre grafeno compacta hexagonal y, (d) P6 sobre grafeno dispo- 


sición entrelazada. 


P3-inter: Complejo mononuclear, donde un átomo central de Co se liga a cuatro molécu- 


las P3. El autoensamblaje, lugar a la interacción entre P3 y formación de una red con com- 
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plejos intercalados (con nanoporos cuadrados). 

P6-Cuadrado: Un átomo de Co central ligado a cuatro moléculas P6. En este complejo 
no hay interacción entre P6. El autoensamblaje da lugar a la formación de una red cuadrada 
con átomos de Co en los vértices. 

P6-inter: Complejo mononuclear, donde un átomo central de Co se liga a cuatro molécu- 
las P6 mediante enlaces Co-NC. El autoensamblaje, da lugar a la interacción entre P6 y 
formación de una red con complejos intercalados (con nanoporos cuadrados). 

Todos los sistemas fueron relajados y se calcularon energías totales. 

Luego de realizar una relajación total de las moléculas en fase gas, la energía total para 
cada configuración fue de -10.7 keV, 5 keV, -16 keV y -30 ke V para P3-Cuadrado, P3-inter, 
P6-Cuadrado y P6-inter, respectivamente. Las estructuras con mínima energía corresponden 
a P3-cuadrado y P6-inter, por tanto, la síntesis experimental es de esperarse. Este resultado 
es acorde a trabajos experimentales sobre el sistema P6-inter, cuya imagen de microscopía 
de barrido de efecto túnel (STM) se muestra en la Fig. (A.2). P3-cuadrado y P6-inter son 
sistemas con número de coordinación 4 para el Co (4:1 ligando/Co), en los dos casos los 
átomos de Co se disponen en una red cuadrada. La diferencia radica en el número de anillos 
feniles de los ligandos y la interacción que surge entre ellos. 

En el caso de P3-cuadrado, las interacciones metal-ligando, son más fuertes que los en- 
laces de hidrógeno y las interacciones tipo van der Waals que se pueden generar entre P3, 
impidiendo la estructura intercalada. 

En P6-intercalado al aumentar el número de anillos feniles, el autoensamblaje se encuen- 
tra dirigido por un lado, por la interacción metal-ligando con coordinación 4 para Co-NC 
y, por el otro, por una interacción lateral no-covalente. Los cálculos DFT muestran que la 
estructura es estable para distancias Co-NC de 1.83 Á. 

Como trabajo futuro se plantea continuar en esta línea, para profundizar en la estructura 
electrónica de éstos materiales de baja dimensionalidad y, determinar propiedades ópticas y 
de transporte electrónico. Específicamente: para los sistemas más estables con Co, agregar la 


superficie de grafeno y calcular la transferencia de carga, energías y distancias de adsorción. 
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application in amino acid sensing 
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Cite this: Phys. Chem. Chem. Phys., 


2019, 21, 597 Sindy J. Rodríguez (3 *? and Eduardo A. Albanesi? 


We modeled a type of field-effect transistor device based on graphene for the recognition of amino 
acids with a potential application in the building of a protein sequencer. The theoretical model used was 
a combination of density functional theory (DFT) with the non-equilibrium Greens function (NEGF) in 
order to describe the coherent transport in molecular devices. First, we studied the physisorption of 
each amino acid on a graphene sheet and we reported the adsorption energy, the adsorption distances, 
the equilibrium configuration and the charge transfer of ten amino acids that can be considered as 
representative of all of the amino acids: histidine (His), alanine (Ala), aspartic acid (Asp), tyrosine (Tyr), 
arginine (Arg), glutamic acid (Glu), glycine (Gly), phenylalanine (Phe), proline (Pro) and lysine (Lys). As a 
result, significant differences were found in the density of states (DOS) after adsorption and there was a 
change in the semi-metallic character of the graphene due to the lysine and arginine interactions. 
Furthermore, we noticed changes in the electrical characteristics of the devices, as the amino acids 
adsorbed onto the surface of the graphene. The curves of current vs. bias voltage (/-V,) display a distinct 
response for each amino acid, ¡.e. the /-V, curves produce a characteristic footprint for each amino 
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acid. We identified a possible rectification mechanism related to the voltage profile asymmetry, where 
the amino acids can control the transport characteristics in the device, ¡.e. Lys and Phe amino acids 
DOI: 10.1039/c8cp05093g physisorbed on graphene act as a molecular diode, where electrons can easily flow in one direction and 
decrease in the other. This may be promising for the prospect of biosensors: graphene could be used as 


rsc.li/pccp an amino acid detector. 


graphene is a candidate for the construction of biosensors at 
the molecular scale. Some advantages of graphene-based bio- 


1 Introduction 


Graphene is a two-dimensional crystal (single layer of carbon 
atoms) that is stable at room temperature and is a gapless 
semiconductor.”” The carbon bonds are sp” hybridized, where 
the in-plane o bonds defined the hexagonal structure and, the 
delocalized tr bonds perpendicular to the planar structure 
are responsible for the electronics properties of graphene.** 
The delocalized electrons give graphene a high electrical con- 
ductivity of 1738 s m *,” and a high charge carrier mobility 
where 15000 cm” V* s * is the mobility measured at room 
temperature for graphene with one or three layers, while a 
single layer achieved 230 000 cm? V”* s * at a temperature near 
0 K.* This endows the graphene with the capacity to form weak 
interactions with other layers of graphene or several other 
substrates. Due to the high electrical sensitivity, the high 
signal-to-noise and the versatility to functionalize the material, 
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sensors are: (i) a high surface-to-volume ratio (2630 m” gy 
which provides high sensitivity to local environmental conditions, 
which is higher than that of graphite (10 m? g *) and CNT 
(1315 m* g *),'** (ii) graphene does not contain metallic 
impurities which interfere with the electrochemistry of the 
material (as occurs with carbon nanotubes CNT), (iii) low 
resistivity with a value of 10? Q cm,'* (iv) low electrical noise 
(and low charge-scattering) and ballistic transport'**” and, 
(v) the ease of interactions with biomolecules by n-1 stacking 
interactions.*” 

Graphene, as a biosensor, provides a variety of detection 
mechanisms, including charge transfer,'” charge scattering 
and field effect. A graphene field effect transistor (GFET) is a 
system with a graphene channel between two electrodes, where 
it is possible regulate the electrical conductivity of the graphene 
upon the application of an external electric field e.g. induced by 
a charged biomolecule (changes in the charge carrier density 
and modulations of the channel conductance are generated 
due to electrical charge redistribution that arises in the 
substrate-molecule interaction process) and further changes 
in current can be observed.'””% In the diversity of theoretical 
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The study of biosensors based on graphene has increased in the last years, the combination of excellent 
electrical properties and low noise makes graphene a material for next generation electronic devices. 
This work discusses the application of a graphene-based biosensor for the detection of amino acids histi- 
dine (His), alanine (Ala), aspartic acid (Asp), and tyrosine (Tyr). First, we present the results of modeling 
from first principles the adsorption of the four amino acids on a graphene sheet, we calculate adsorp- 
tion energy, substrate-adsorbate distance, equilibrium geometrical configurations (upon relaxation) and 





o a densities of states (DOS) for each biomolecule adsorbed. Furthermore, in order to evaluate the effects of 
Biosensof amino acid adsorption on the electronic transport of graphene, we modeled a device using first-principles 
Amino acid calculations with a combination of Density Functional Theory (DFT) and Nonequilibrium Greens Func- 
Adsorption tions (NEGF). We provide with a detailed discussion in terms of transmission, current-voltage curves, 


and charge transfer. We found evidence of differences in the electronic transport through the graphene 
sheet due to amino acid adsorption, reinforcing the possibility of graphene-based sensors for amino acid 


sequencing of proteins. 


O 2017 Elsevier B.V. All rights reserved. 





1. Introduction 


The analysis and quantification of biomolecules is crucial in clin- 
ical diagnosis and treatment, for this reason in the last years, the 
construction of biosensors with biomedical application has gained 
great importance [1-3]. For building biosensors, it is needed to 
explore materials with high biocompatibility, sensitivity, selec- 
tivity with fast response time, and feasible nanoscale fabrication 
procedures. Graphene, a single layer of carbon atoms, has exhibited 
superior physical and chemical properties than other 3D materials, 
positionning itas a strong candidate for the construction of biosen- 
sors [4,5]. This material is characterized as a semi-metal or zero 
gap semiconductor. As for its electrical properties, it has shown (1) 
a remarkably high electron mobility at room temperature—with 
experimentally reported values in excess of 15,000 cm? V-1 s71 
[6]—, (11) low resistivity (1078 (2 cm), (iii) low Johnson noise, which 
along with its high electron mobility allow it to be utilized as 
the channel in a field effect transistor (FET), (iv) high surface area 


* Corresponding author. 
E-mail address: sindy.rodriguezQifis.santafe-conicet.gov.ar (S.J. Rodríguez). 
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0169-4332/0 2017 Elsevier B.V. All rights reserved. 


2620 m?/g for both sides of graphene [7], and 1310 m?/g for one- 
side (e.g., supported on a substrate). 

Theoretical and experimental advances in structures of 
graphene-based nanomaterials reported changes in electronic 
transport properties of a graphene sheet, due to interactions by 
covalent or non-covalent forces between graphene and sever- 
als organic molecules [8-12]. Viswanathan et al. [13] described 
an approach for the development of a graphene-based biosen- 
sor platform using glucose as an example of target molecule. 
The presence of external molecules can vary its conductivity and 
this variation can either be monitored using a simple chemiresis- 
tor or by a transistor based sensor. Ohno et al. [14] investigated 
graphene field-effect transistors (GFETs) for electrical detection of 
pH and protein adsorptions, the GFETs thus acted as highly sen- 
sitive electrical sensors for detecting biomolecule concentrations. 
Furthermore, smaller molecules have been sensed: dopamine [15], 
and nucleotides in a DNA chain, among others. Zou et al. [16] 
reported a DNA sensor based on graphene, the current signals of 
the four bases guanine (G), adenine (A), thymine (T) and cyto- 
sine (C), were separated efficiently. Zhen et al. [17] developed a 
novel FET nanobiosensor based on a chemical vapor deposition 
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Abstract The effects of histidine and its imidazole ring 
adsorption on the electronic transport properties of graphene 
were investigated by first-principles calculations within a 
combination of density functional theory and non-equilib- 
rium Greens functions. Firstly, we report adsorption energies, 
adsorption distances, and equilibrium geometrical configu- 
rations with no bias voltage applied. Secondly, we model 
a device for the transport properties study: a central scat- 
tering region consisting of a finite graphene sheet with the 
adsorbed molecule sandwiched between semi-infinite source 
(left) and drain (right) graphene electrode regions. The elec- 
tronic density, electrical current, and electronic transmission 
were calculated as a function of an applied bias voltage. 
Studying the adsorption of the two systems, 1.e., the his- 
tidine and its imidazole ring, allowed us to evaluate the 
importance of including the carboxyl (-COOH) and amine 
(NH) groups. We found that the histidine and the imida- 
zole ring affects differently the electronic transport through 
the graphene sheet, posing the possibility of graphene-based 
sensors with an interesting sensibility and specificity. 


Keywords Adsorption amino acid - Graphene - NEGF 
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1 Introduction 


Graphene is a semimetal with zero energy bandgap, and the 
carbon atoms are arranged in a hexagonal lattice in a single- 
atom-thick planar sheet. In comparison with other 3D mate- 
rials, graphene has unique physical and chemical properties: 
high carrier mobility (can exceed 15,000 cm? V-! 571 
even under ambient conditions), extremely high surface-to- 
volume ratio, flexibility, transparency, and biocompatibility, 
among others [1]. These features make it a good candidate 
for a new generation of electronic devices and bionic tech- 
nology. In particular, the construction of nanobiosensors has 
become increasingly relevant because many biological pro- 
cesses and pharmacological responses may be monitored 
1f biomolecules are detected and quantified [2]. Sensors 
for biological application must be biocompatible, sensitive, 
selective, and also with excellent electrical properties (such 
as a low intrinsic electronic noise): graphene could play an 
important role in this area. For instance, Baraket et al. [3] 
showed experimentally that aminated graphene can work as 
a biologically active field-effect transistor for DNA detec- 
tion. Also, several theoretical and experimental studies have 
proposed the use of graphene in the construction of biosen- 
sors [2-7], e.g., graphene to detect glucose [8], proteins [9], 
DNA [10], or even bacteria [11,12]. 

Several biomolecules tend to bind to graphene through 
non-covalent bonds, such as m—1T stacking interaction. These 
interactions are relevant in the design and fabrication of 
nanodevices, because subtle changes in the electronic charac- 
teristics of the 71 systems with graphene can lead to relevant 
defects in the nanosystem structure and properties [5]. The- 
oretical studies based on ab initio calculations have been 
instrumental in understanding the nature of —T interactions. 
Song et al. [13] reported that oxygen in nucleobases adsorbed 
on graphene with 1—7x stacking interaction alters the elec- 
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Abstract. Previous studies have demonstrated how the interactions between biomolecules and 
graphene play a crucial role in the characterization and functionalization of biosensors. In this 
paper we present a theoretical study of the adsorption of histidine on graphene using density 
functional theory (DFT). In order to evaluate the relevance of including the carboxyl (-COOH) 
and amino (-NH2) groups in the calculations, we considered i) the histidine complete (i.e., with 
its carboxyl and its amino groups included), and ii) the histidine's imidazole ring alone. We 
calculated the density of states for the two systems before and after adsorption. Furthermore, we 
compared the results of three approximations of the exchange and correlation interactions: local 
density (LDA), the generalized gradients by Perdew, Burke and Ernzerhof (GGA-PBE), and 
one including van der Waals forces (DFT-D2). We found that the adsorption energy calculated 
by DF'T-D2 is higher than the other two: Esas prr-D2 >Eaas-LDaA >Eads-aGGA. We report 
the existence of charge transfer from graphene to the molecule when the adsorption occurs; this 
charge transfer turns up to be greater for the complete histidine than for the imidazole ring 
alone. Our results revealed that including the carboxyl and amino groups generates a shift in 
the states of imidazole ring towards Ep. 


1. Introduction 

Graphene is a two-dimensional crystal with special chemical and physical properties [1, 2]. 
Consequently, since its discovery, it has been used in fields such as electronics, biomedicine, 
biotechnology and pharmacology, among others [3, 4, 5]. For instance, Wang et al. [6], in the 
area of drug delivery devices and molecular transport, reported that graphene could be used 
as a good vehicle of DNA aptamers into living cells, because it provides effective protection 
of the oligonucleotides from enzymatic cleavage. Furthermore, Hess et al. [7, 8] developed 
graphene transistor arrays for recording cellular action potentials, which have applications in 
neural prosthesis. The study suggests the simple integration of graphene electronics with flexible 
technologies (an important aspect of biomedical implants design), because they prevent severe 
tissue damage [7]. 

In the area of molecular biology and nanomaterials, different studies have been developed on 
the adsorption of biological molecules on graphene sheets and carbon nanotubes [9, 10, 11, 12]. 
Graphene sensors can provide high sensitivity, selectivity and low intrinsic electronic noise, which 
is a great advantage over other materials [3]. Chemical sensors of carbon nanotubes for small 


Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution 
BY of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOL. 
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Q2 Abstract 


The construction of sensors plays an important role in electronics and bioelectronics because of 
their wide range of novel applications for biomedical research, food quality control, and environ- 
mental monitoring. The sensors for biological application demand materials with special prop- 
erties, such as: sensitivity, selectivity, biocompatibility, high electronic mobility, low electronic 
noise, and chemical functionality. Due to its unique physical and chemical properties, graphene 
has emerged as a suitable candidate for making sensors. "The electronic properties of graphene 
are extremely sensitive to environmental perturbations such as electronic doping and molecu- 
lar adsorption. Several biomolecules interact with graphene through noncovalent bonds, and 
these interactions can modify the electron density and conductivity in the graphene allowing the 
detection of the molecules. In the fabrication of electronic components like field-effect transistors 
(FE Ts) with the adoption of thin layers of graphene, which consist of two terminals, the source 
and drain, and a gate that controls the resistance of the device, the highly mobile electrons at or 
near its surface are extremely sensitive to local charge changes. As a result, the molecules acquire 
charge when adsorbed, and their binding to a graphene-based gate will disrupt the flow of the elec- 
trons. Studying and quantifying the electronic eftects due to interaction between substrates and 
adsorbates can provide tools for the construction of nanographene-based sensing devices. But, 
what are the phenomena that underlie the charge transfer between graphene and biomolecules? 
The different effects generated in the adsorption process of biomolecules and the modification and 
response of the electrical properties of the graphene, are presented in this review. "The chapter has 
four sections: i) introduction, where the generalities and basic notions of the electronic properties 
for graphene are presented, together with a brief fundamental of electronic transport based on 
nonequilibrium Greens functions, and the experimental approaches to current-voltage curves in 
FET devices; ii) computational modeling of adsorption and devices; iii) experimental realization 
of devices; and iv) conclusions and final remarks. 
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Abstract— The possibility of utilizing graphene as a biosensor 
has received a considerable attention in recent years. Here, we 
present results on the molecular dynamics of the adsorption of 
the aminoacid histidine on graphene, based on ab initio calcu- 
lations within a pseudopotentials approach. Taking into consid- 
eration the amine and carboxylic groups of the aminoacid (not 
only the imidazole ring), we calculate the adsorption energy and 
final mean distance to the graphene sheet, along with the slight 
deformation of the graphene. Furthermore, we provide with a 
detailed discussion on two ways of calculating (and presenting) 
the adsorption curve (interaction energy vs. mean distance), de- 
pending on how relaxation of atomic positions is included in the 
calculations. 


Keywords— graphene-based biosensors, hisidine aminoacid, 
Density Functional Theory. 


Il. INTRODUCTION 


Graphene is one very particular allotropic form of car- 
bon (the others being amorphous carbon, graphite, and di- 
amond), since it consists of a single atomic layer, i.e., 1t is 
a two-dimensional material. Even though this material had 
been theoretically known and studied for decades, it has re- 
ceived an enormous amount of attention since the works of 
A. Geim [1] and K. Novoselov [2], who received the Nobel 
Prize in 2010 for their experimental achievements regarding 
the material. 

The possibility of using graphene as a chemical sensor has 
been studied for several years now (e.g., Refs. [3, 4]) and it 
1s based on the fact that graphene is intrinsically a low-noise 
material [5]. While the adsorption of histidine on graphene, 
along with the other three aromatic aminoacids (phenylala- 
nine, tyrosine, and tryptophan), has already been studied by 
ab initio calculations [6], in that work only the imidazole ring 
was taken into consideration; while that is, indeed, the dis- 
tinctive residue of histidine, the molecule in its physiologi- 
cal form does also have an amino and a carboxylic group. 
In fact, the latter has been reported to bear a critical effect 
on the electronic properties of graphene near the Fermi en- 
ergy [7]. Hence, we decided to undertake the study of the 
adsorption process considering the complete molecule (in its 


O Springer International Publishing Switzerland 2015 


non-protonated state) and focusing on the energetics and the 
geometrical changes during the process of adsorption. 


II. COMPUTATIONAL DETAILS 


Figure 1 shows a molecule of histidine at a mean distance 
h from a sheet of graphene (periodic conditions make the 
graphene sheet infinite). The latteral dimensions chosen for 
the supercell (12.78 x12.29 x25 Á) guarantee the decoupling 
of adjacent occurrences of the molecule, while the height of 
25 Á allows us to locate the molecule far enough from the 
graphene, so a reference energy (E"*/) can be determined for 
the latter calculation of the interaction energy (attraction 1f far 
enough, and rejection if too close): E" (h) = E'%(h) — Ef, 
where E"*f = E'"(h = 12.5 Á). The adsorption energy is that 
of the interaction at the global minimum of the curve E”" (h). 





Fig. 1: The histidine molecule on graphene. 


Al the calculations were performed within a pseudopo- 
tentials approach to the Density Functional Theory, utilizing 
the OpenMX [8] software suite and adopting a Generalized 
Gradient Approximation (PBE-GGA [9]) for the exchange- 
correlation potential. A cut-off energy of 180 Ry and a k- 
mesh of 5x3 x1 were used for the self-consistency, and for 
the calculation of the densities of states (DOS) the grid was 
refined up to 50 x 50 x 1. For the relaxation, the convergence 
criterion was of 0.0004 Ha/Bohr. 


140 
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Resumen 


El grafeno es una forma alotrópica del carbono, sintetizado por primera vez por Geim y 
Novoselov en 2004, lo que les valió el Premio Nobel en el año 2010. Este material conforma una 
red cristalina bidimensional de un sólo átomo de espesor, siendo la disposición espacial de sus 
átomos semejante a un panal de abejas. Presenta extraordinarias propiedades electrónicas y 
fisicoquímicas como las de tener portadores de carga con masa efectiva cero con altísima 
movilidad, coeficiente térmico negativo, ser flexible, transparente, liviano y de extrema dureza, 
entre otras. 


Las múltiples propiedades del grafeno le han abierto una amplia gama de aplicaciones en 
áreas como la nanotecnología, biomedicina, farmacología, bioelectrónica, entre otras. 
Particularmente, en la construcción de biosensores, el grafeno es un material que ofrece 
muchas ventajas, ya que, es sensible, selectivo, de bajo ruido electrónico intrínseco y además, 
biocompatible. Trabajos experimentales y teóricos dan cuenta de la facilidad del grafteno para 
adsorber moléculas biológicas, hecho que permite estudiar qué propiedades específicas 
modifica el material luego de la adsorción, información útil para la eventual funcionalización y 
caracterización de un biosensor basado en grafeno. Trabajos en esta dirección, con 
nuecleobases de ADN han sido explorados por diversos investigadores. 


En este trabajo, se modelizó mediante técnicas basadas en la teoría de la funcional densidad 
(DFT) la adsorción de un aminoácido sobre grafeno, específicamente, histidina. Se realizaron 
gráficas de densidades de estados y densidades de carga para el sistema antes y después de 
la adsorción, con el fin de identificar que átomos de la molécula y del grafeno tienen mayor 
participación en la redistribución de estados y carga. En los resultados reportados para la 
energía y distancia de adsorción se compara, por un lado, diferencias entre tres esquemas 
para la interacción de intercambio y correlación entre electrones: de densidad local (LDA]); de 
gradientes generalizados (GGA-PBE): y de van der Waals (DFT-D2). Por otro lado se estudió la 
relevancia de la inclusión en los cálculos, de los grupos carboxilo y amino comunes a todos los 
aminoácidos. A futuro se evaluarán los efectos de la adsorción en el transporte electrónico en 
el grafeno. 
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Resumen 


Actualmente, el autoensamblaje de moléculas ha mostrado ser una estrategia efectiva para manipular la propiedades electrónicas de metales. 
Las moléculas autoorganizadas en una red porosa, funcionan como una barrera de potencial que es capaz de confinar los electrones libres de 
la superficie de los metales, cambiando así, su estructura de banda. Este enfoque, ha inspirado el estudio de moléculas orgánicas en grafeno, 
como una posible ruta para modificar sus propiedades eléctricas. Sin embargo, a pesar de la importancia de las redes porosas en grafeno, aún 
hace falta un claro entendimiento del autoensamblaje de estas moléculas en grafeno. En este trabajo presentamos los resultados de la 
modelización del autoensamblado de moléculas orgánicas de Di-carbonitrilo-polifenil (NC-Ph,-CN, donde n=3 y 6) y cobalto (Co) sobre una 
hoja de grafeno. Los resultados teóricos desde la teoría de la funcional densidad (DFT), dan cuenta de la diferencia de periodicidad y 
competencia entre las interacciones no covalentes que dirigen el autoensamblado de redes, al modificar las longitudes de las moléculas 


orgánicas (NC-Ph,-CN y NC-Ph¿-CN con longitudes de 16.8 A y 29.6 A respectivamente). Todas las redes modeladas tienen un número de 
coordinación 4, con centros de cobalto en los nodos. Se reportan los resultados de energía de enlace, transferencia de carga y distancias de 


adsorción. Los resultados teóricos que modelan las configuraciones adoptadas, son acordes a los resultados experimentales y resultan un 
avance considerable en nuestro entendimiento del autoensamblaje molecular en grafeno. 


Palabras clave: Grafeno, autoensamblaje, polifenil. 
Keywords: Graphene, self-assembly, polyphenol. 


1. Introducción 

El autoensamblaje molecular sobre superficies planas es un 
campo que ha despertado un creciente interés en los últimos 
años. Recientemente, en el 2016 Jean-Pierre Sauvage, J. Fraser 
Stoddart y Bernard L. Feringa, reciben el premio nobel de 
química, gracias al diseño y la síntesis de máquinas 
moleculares. La construcción de máquinas moleculares y de 
redes porosas por autoensamblaje, tienen aplicaciones en el 
desarrollo de nuevos materiales con múltiples propiedades 
Ópticas y eléctricas, sensores y nuevos sistemas de 
almacenamiento de energía. 

El autoensamblaje molecular es el proceso de asociación 
espontánea de moléculas para la formación de estructuras más 
complejas como supramoléculas, estructuras porosas, entre 
otras [1]. La formación de las redes se encuentra dirigido por 
las interacciones molécula-molécula y, molécula-substrato. 
Cuando las moléculas son depositadas sobre superficies poco 
reactivas, como grafeno, la estabilidad, dimensión y geometría 
de los poros, depende principalmente de la interacción entre 


moléculas de la capa, ya que son relativamente fuertes y 
altamente  direccionadas, abarcando interacciones no 
covalentes, metal-ligando y enlaces de hidrógeno [2]. 

En este trabajo se llevó a cabo un modelado computacional 
desde la teoría del funcional de la densidad (DFT), de la 
formación de redes autoensambladas con ligandos de 
Di-carbonitrilo-polifenil NC-Ph,-CN (P6), NC-Ph,-CN (P3) y 
cobalto (Co) sobre un substrato de grafeno. Se estudiaron: 1) 
cuatro configuraciones moleculares de ligandos y Co en fase 
gas, se discute en términos de energía total y densidad de carga 
la posibilidad de su síntesis experimental y, ii) dos 
configuraciones con substrato de  grafeno, calculando 
transferencia de carga, energías y distancias de adsorción. Los 
resultados teóricos permiten comprender, por un lado, el rol 
que juega el número de anillos fenilos en la geometría de los 
poros, y por el otro, las interacciones que dirigen el 
autoensamblado sobre grafeno. 

A futuro, se espera funcionalizar el grafeno, de tal forma que, 
las redes porosas autoensambladas sirvan como máscaras para 


ETERNA] 


97 


Bibliografía 


[1] Novoselov, K. and Geim, A. and Morozov, $. and Jiang, D. and Zhang, Y. and Dubonos, 
S. and Grigorieva, I. and Firsov, A. Two-dimensional gas of massless Dirac fermions 
in graphene. Nature, 438(7065):197—200, 20053. 


[2] Toan, T.T and Mulchandani, A. Carbon nanotubes and graphene nano field-effect 
transistor-based biosensors. T.-T. Tran, 76(1):222-232, 2016. [xvii] [9] [10] 


[3] Mohanty, N. and Berry, V. Graphene-based single-bacterium resolution biodevice and 
DNA transistor: Interfacing graphene derivatives with nanoscale and microscale bio- 
components. Nano Letters, 8(12):4469-4476, 2008. 


[4] Hess, L. and Jansen, M. and Maybeck, V. and Hauf, M. and Seifert, M. and Stutzmann, 
M. and Sharp, L. and Offenhusser, A. and Garrido, J. Graphene transistor arrays for re- 


cording action potentials from electrogenic cells. Advanced Materials, 23:3043-3049, 


2011. xvi] [11] 


Singh, D. and Nordstrom, L. Planewaves, Pseudopotentials, and the LAPW Method. 


Springer, 2006. 


Singla, P. and Ryaz, M. and Singhal, S. and Goel, N. Theoretical study of adsorption of 


[S 


| 


[6 


La) 


amino acids on graphene and BN sheet in gas and aqueous phase with empirical DFT 
dispersion correction. Phys. Chem. Chem. Phys., 18(7):5597-5604, 2016. 


53/154)[55 


[7 


La) 


Tavassoli, L. and Darvish, M. and Jahanshahi, M. Trends of amino acid adsorption onto 


graphene and graphene oxide surfaces: a dispersion corrected DFT study. RSC Adv., 


5:92843-92857, 2015. 


99 


BIBLIOGRAFÍA 


[8] Jiang, D. and Chen, Z. Graphene Chemistry: theoretical perspective. John Wiley and 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


Sons, Ltd, United Kingdom, 2013. 16] 


Georgakilas, V. and Otyepka, M. and Bourlinos, A. and Chandra, V. and Kim, N. and 


Kemp, K. and Hobza, P. and Zboril, R. Funcionalization of graphene: Covalent and 


Non-Covalent Approaches, Derivatives and Applications. Chem. Rev., 112:6156-6214, 


2012. lxx1)/41[8] 


Hess, L.H. and Se1fert, M. and Garrido, J.A. Graphene transistors for bioelectronics. 


Proc. IEEE, 101:1780-1792, 2013. 


Ferrari, A.C et al. Science and technology roadmap for graphene, related two- 
dimensional crystals, and hybrid systems. Nanoscale, 7:4598, 2015. 


Novoselov, K. and Geim, A. and Morozov, S. and Jiang, D. and Zhang, Y. and Dubonos, 
S. and Grigorieva, I. and Firsov, A. Electric Field Effect in Atomically Thin Carbon 
Films. science, 306:666-668, 2004. 


Tuan, D. Charge and Spin Transport in Disordered Graphene-Based Materials. Sprin- 
ger, 2016. 


Mayorov, A.S. and Gorbachev, R.V. and Morozov, S.V. and Britnell, L. and Jalil, R. and 
Ponomarenko, L.A. and Blake, P. and Novoselov, K.S. and Watanabe, K. and Taniguchi, 
T. and Geim, A.K. Micrometer-Scale Ballistic Transport in Encapsulated Graphene at 
Room Temperature. Nano. Lett, 11:2396-2399, 2011. 


Tse, W.K. and Hwang, E.H. and Das Sarma, S. Ballistic Hot Electron Transport in 
Graphene. Appl. Phys. Lett, 93:023128, 2008. 


Weiss, N. and Zhou, H. and Liao, L. and Liu, Y. and Jiang, S. and Huang, Y. and Duan, 
X. Graphene: An emerging electronic material. Advanced Materials, 24(43):5782- 
5825, 2012. 


Zhu, Y. and Mural, S. and Cai, W. and Li, X. and Suk, J.W. and Potts, J.R. and Ruoff, 
R.S. . Graphene and graphene oxide: Synthesis, properties, and applications. Ad». 
Mater, 22:3906-3924, 2010. 


100 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


BIBLIOGRAFÍA 


Liu, Y. and Dong, X. and Chen, P. Biological and chemical sensors based on graphene 
materials. Chem Soc Rev, 41(6):2283-2307, 2012. 


Bolotin, K.I. and Sikes, K.J. and Jiang, Z. and Klima, M. and Fudenberg, G. and Hone, 
J. and Kim, P. and Stormer, H.L. Ultrahigh electron mobility in suspended graphene. 
Solid State Communications, 146(9-10):351-355, 2008. 


Amin, R.K. and Bid, A. Graphene as a sensor. CURRENT SCIENCE, 107(c):430-433, 
2014. 


Saha, B. and Bhattacharyya, P. Adsorption of amino acids on boron and/or nitrogen 
doped functionalized graphene: A Density Functional Study. Computational and Theo- 
retical Chemistry, 1086:45-51, 2016. 


Mallakpoura, S. and Abdolmalekia, A. and Borandeha, S. Covalently functionalized 
graphene sheets with biocompatible natural amino acids. Applied Surface Science, 
307(307):533-542, 2014. 


Ghosh, A. and Rao, K. and George, S.J. and Rao, C.N. Noncovalent Functionaliza- 
tion, Exfoliation, and Solubilization of Graphene in Water by Employing a Fluorescent 
Coronene Carboxylate. Chem. Eur. J, 16:2700 -2704, 2010. [6] 


Cheng, H. and Shiue. R. and Tsai1, C. and Wang, W.H. and Chen, Y. High-Quality 
Graphene p-n Junctions via Resist-free Fabrication and Solution-Based Noncovalent 
Functionalization. AcsNANO, 5:2051-2059, 2011. [6] 


Gowtham, $. and Scheicher, R.H. and Ahuja, R. and Pandey, R. and Karna, S.P. Phy- 
sisorption of nucleases on graphene: Density-functional calculations. Physical Review 
B, 76:033401, 2007. 


Lee, E.C. Effects of DNA nucleotide adsorption on the conductance of graphene nano- 
ribbons from first principles. Applied Physics Letters, 100:153117, 2012. 


Le, D. and Kara, A. and Schróder, E. and Hyldgaard, P. and Rahman, T. Physisorption 
of nucleobases on graphene: a comparative van der Waals study. Journal of Physics: 
Condensed Matter, 24(42):424210, 2012. 


101 


BIBLIOGRAFÍA 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[39] 


[35] 


[36] 


[37] 


Zhiani, R. Adsorption of various types of amino acids on the graphene and boron- 
nitride nano-sheet, a DFT-D3 study. Applied Surface Science, 409:35-—44, 2017. 


Rodríguez, S.J. and Makinistian, L. and Albanesi, E.A. Graphene for amino acid 
biosensing: Theoretical study of the electronic transport. Applied Surface Science, 
419:540-545, 2017. 


Fernandez, A.C.R. and Castellani, N.J. Noncovalent Interactions between Dopamine 
and Regular and Defective Graphene. ChemPhysChem, 2017. 


Ortiz-Medina, J. and López-Urías, F. and Terrones, H. and Rodríguez-Macías, F.J. and 
Endo, M. and Terrones, M. Differential response of doped/defective graphene and 
dopamine to electric fields: A density functional theory study. Journal of Physical 
Chemistry C, 119(24):13972-13978, 2015. 


Liu, S.J. and Wen, Q. and Tang, L.J. and Jiang, J.H. Phospholipid-Graphene Nanoas- 
sembly as a Fluorescence Biosensor for Sensitive Detection of Phospholipase D Acti- 
vity. Anal. Chem, 84:5944-5950, 2012. 


Hu, Y. et al. Biocompatible Graphene for Bivanalytical Applications. SpringerBriefs 
in Molecular Science, 2015. 


Gan, X.R. and Zhao, HM. A Review: Nanomaterials Applied in Graphene-Based 
Electrochemical Biosensors. Sensors and Materials, 2702):191-215, 2015. 


Zhan, B. and Chen, L. and Jun, Y. and Gareth, J. and We1, H. and Xiaochen, D. Graphe- 
ne field-effect transistor and its application for electronic sensing. Small, 10(20):4042— 
4065, 2014. 


Yang, W. and Ratinac, K. and Ringer, S. and Thordarson, P. and Gooding, J. and Braet, 
F. Carbon nanomaterials in biosensors: Should you use nanotubes or graphene. Ange- 
wandte Chemie - International Edition, 49(12):2114-2138, 2010. 


Viswanathan, S. and Narayanan, T. and Aran, K. and Fink, K. and Paredes, J. and Aja- 
yan, P. and Filipek, S. and Miszta, P. and Tekin, C. and Inc1, F. and Demirci, U. and Li, 
P. and Bolotin, K. and Liepmanmn, D. and Renugopalakrishanan, V. Graphene-protein 
field effect biosensors: Glucose sensing. Materials Today, 18(9):513-522, 2015. 


102 


[38] 


[39] 


[40] 


[41] 


[42] 


[43] 


[44] 


[45] 


[46] 


[47] 


[48] 


[49] 


BIBLIOGRAFÍA 


Yan, FE. and Zhang, M. and Li, J. Solution-Gated Graphene Transistors for Chemical 
and Biological Sensors. Advanced Healthcare Materials, 3(3):313-331, 2014. [8] [9] 


Reddy, D. and Register, L. and Carpenter, G. and Banerje, S.K. Graphene field-effect 
transistors. Journal of Physics D: Applied Physics, 45(1):019501, 2012. [9] 


Green, N.S. and Norton, M. Interactions of DNA with graphene and sensing appli- 
cations of graphene field-effect transistor devices: A review. Analytica Chimica Acta, 
853(1):127-142, 2015. [9] 


He, Q. and Wu, $. and Yin, Z. and Zhang, H. Graphene-based electronic sensors. Chem. 
Sci., 3(1):1764-1766, 2012. 


Bendali, A. and Hess, L. and Se1fert, M. and Forster, V. and Stephan, A. and Garrido, J. 
and Picaud, S. Purified neurons can survive on peptide-free graphene layers. Advanced 
Healthcare Materials, 2:929-933, 2013. 


Schedin, F. and Geim, A.K. and Morozov, S.V. and Hill, E.W. and Blake, P. and Kats- 
nelson, M.I. and Novoselov, K.S. . Detection of individual gas molecules adsorbed on 
eraphene. nature materials, 6:652-635, 2007. 


Krupski, J. PPS (b)., 157:199—207, 1990. 


Patil, A.B. and Fernandes, F.B. and Bueno, P.R. and Davis, J. Electrochemical sensing 
of L-histidine based on structure-switching DNA zymes and gold nanoparticle-graphene 
nanosheet composites. Bioanalysis, 71:723—742, 2015. 


Song, B. and Cuniberti, G. and Sanvito, S. and Fang, H. Nucleobase adsorbed at grap- 
hene devices: Enhance bio-sensorics. Appl. Phys. Lett, 100:063101, 2012. 


Springborg, M. Density-Functional Methods in Chemistry and Materials Science. Wi- 


ley, New York, 1997. 


Kohn, W. Nobel lecture: Electronic structure of matter?ave functions and density fun- 
ctionals. Rev Mod Phys., 71:1253-1266, 1999. 


P. Hohenberg and W. Kohn. Inhomogeneous electron gas. Phys Rev., 136:B864—B871, 
1964. 


103 


BIBLIOGRAFÍA 


[50] Kohanoff, J. Electronic Structure Calculations For Solids and Molecules: Theory and 
Computational Methods. Cambridge University Press, 2006. 


[51] Woods, L.M. and Badescu, C. and Reinecke, T.L. Adsorption of simple benzene der1- 
vatives on carbon nanotubes. Physical Review B, '75:153415-1535425(9), 2007. 


[52] Ortmann, F. and Schmidt, W.G. and Bechstedt, FE Atracted by long-range electron 
correlation: Adenine on graphite. Physical Review Letters, 95:186101-186104, 2005. 


[53] Grime, S. and Antony, J. and Ehrlich, S. and Krieg, H. A consistent and accurate 
ab initio parametrization of density functional dispersion correction (dft-d) for the 94 
elements h-pu. J. Chem. Phys, 132:154104, 2010. 


[54] Grime, $. and Ehrlich, S. and Goerigk, L. Effect of the damping function in dispersion 
corrected density functional theory. J. Comput. Chem, 32:1456, 2011. 


[55] Becke, A.  Density-functional thermochemistry. v. systematic optimization of 
exchange-correlation functionals. Y Chem Phys, 107:8544(20), 1997. 


[56] Pomorski, P. and Odbadrakh, K. and Sagui, C. and Roland, C. Nonequilibrium's fun- 
ction modeling of the quantum transport of molecular electronic devices. Theor. Com- 
put. Chem., 17:187-204, 2007. 


[57] Keldysh, L.V. Diagram Technique for Nonequilibrium Processes. Sov. Phys. JETP, 
20:1018, 1965. 


[58] Ozaki, T and Nishio, K. and Kino, H. Efficient implementation of the nonequilibrium 
green function method for electronic transport. Phys Rev B, 81:035116, 2010. 


[59] Pastawsk1, H. and Medina, E. Tight binding methods in quantum transport through 
molecules and small devices: form the coherent to the decoherent description. efficient 
implementation of the nonequilibrium green function method for electronic transport. 
Revista mexicana de Fisica, 47:1-23, 2001. 


[60] Ozaki, T. Numerical atomic basis orbitals from h to kr. Phys. Rev. B., 67:155108, 2003. 
201 


104 


[61] 


[62] 


[63] 


[64] 


[65] 


[66] 


[67] 


[68] 


[69] 


[70] 


BIBLIOGRAFÍA 


Monkhorst, J.H and Pack, J.D. Special points for brillonin-zone integrations. Conden- 
sed Matter and Materials Physics, 13:5188-5192, 1976. 


Areshkin, D.A. and Nikolie, B.K. . Electron density and transport in top-gated graphene 
nanoribbon devices: First-principles green function algorithms for systems containing 
a large number of atoms. Phys. Rev. B, 81:155450, 2010. 


Landauer, R. Spatial variation of currents and fields due to localized scatterers in me- 
tallic conduction. IBM J. Res. Dev., 223:1018, 1957. 


Meir, Y. and Wingreen, N.S. Landauer formula for the current through an interacting 
electron region. Phys. Rev. Lett., 68:2512-2515, 1992. 


Remko, M. and Fitz, D. and Rode, B.M. Effect of Metal lons (Li+, Na+, K+, Mg2+, 
Ca2+, Ni2+, Cu2+, and Zn2+) and Water Coordination on the Structure and Properties 
of L-Arginine and Zwitterionic L-Arginine. J. Phys. Chem. A, 112:7652—7661, 2008. 


Chapo, C.J and Paul. J.B. and Provencal, K. and Saykally R.J. Is Arginine Zwitterionic 
or Neutral in the Gas Phase? Results from IR Cavity Ringdown Spectroscopy. J. Am. 
Chem. Soc., 120:12956-12957, 1998. 


Rodríguez, S.J. and Makinistian, L. and Albanesi, E. Molecular dynamics simulation 
of the adsorption of histidine on graphene. IFMBE Proc. Springer Series, 49:140-143, 


2015. 


Rodríguez, S.J. and Makinistian, L. and Albanesi, E. Theoretical study of the adsor- 
ption of histidine amino acid on graphene. Journal of Physics, 705:17-19, 2016. 


Rodríguez, S.J. and Makinistian, L. and Albanesi, E.A. Computational study of trans- 
port properties of graphene upon adsorption of an amino acid: importance of including 
-NH) and -COOH groups. J. Comput Electron, 16:127-132, 2017. 


Goerigk, L. A Comprehensive Overview of the DFT-D3 London-Dispersion Correction. 
John Fedor, 1 edition, 2017. 


105 


BIBLIOGRAFÍA 


[71] Reimers, J. and Li, M. and Wan, D. and Gould, T. and Ford, M. . Surface Adsorption. 


[72] 


[73] 


ña 


[75] 


[76] 


[77] 


[78] 


[79] 


[80] 


John Fedor, 1 edition, 2017. 


Ozaki, T. and Kino, H. Ground state of the electron gas by a stochastic method. Phys. 
Rev. B., 69:195113, 2004. 


Berman, H. and Westbrook, J. and Feng, Z. and Gilliland, G. and Bhat, T. and Weissig, 
H. and Shindyalov, I. and Bourne, P. The Protein Data Bank. p, Nucl Acids Res. 28. 
235-242., 2000. 


Tao, N.J. and Shi, Z. Monolayer guanine and adenine on graphite in nacl solution: A 
comparative. J. Phys. Chem., 98:1464-1471, 1994. 


Rangel, N. and Seminario, J. Vibronics and plasmonics based graphene sensors. J. 
Chem. Phys., 132:125102, 2010. 


Liu, E. and Li, X.F. and Yan, Q. and Li, Q. and Zhang, X. Deng, M. and Qiua, Q. and 
Luo, Y. Uniform and perfectly linear current—voltage characteristics of nitrogen-doped 


armchair graphene nanoribbons for nanowires . Phys. Chem. Chem. Phys., 19:44-48, 


2017. 


Zhang, Y. and Zhou, K. and Xie K. and Zeng J. and Zhang, H. and Peng, Y. Tuning the 
electronic structure and transport properties of graphene by noncovalent functionaliza- 


tion: effects of organic donor, acceptor and metal atoms. Nanotechnology, 21:065201 


(Tpp), 2010. 


Kvashnin, A.G. and Kvashnin D.G. and Kvashnina O.P. and Chernozatonskii, 


L.A. Transport investigation of branched graphene nanoflakes. Nanotechnology, 
26:385705(6pp), 2015. 


Sun, J. and Lin, N. and Li, Z. and Ren, H. and Tanga C. and Zhao, X. Electronic and 
transport properties of graphene with grain boundaries. RSC Ad»v., 6:1090-1097, 2016. 


Murali, R. and Yang, Y. and Brenner, K. and Beck, T. and Meindi J.D. Breakdown 
current density of graphene nanoribbons. Appl. Phys. Lett., 94:243114, 2009. 


106 


[81] 


[82] 


[83] 


[84] 


[85] 


[86] 


[87] 


[88] 


[89] 


BIBLIOGRAFÍA 


Castro, N. and Peres, N.M. and Novoselov, K.S. and Geim, A.K. The electronic pro- 
perties of graphene. Rev. Mod. Phys., 81:109-162, 2009. 


Krebs, E. and Grabill, L. and Riemann. Amino acid nanopatterning on graphite. Surface 
Science, 1(1):16604-16611, 2018. 


Pandit, S. and De, M. Interaction of Amino Acids and Graphene Interaction of Amino 
Acids and Graphene Oxide : Trends in Thermodynamic Properties. J. Phys. Chem. C, 
121:600—608, 2016. 


Zhang, M. and Yin, B. and Wang, X. and Ye, B. Interaction of peptides with graphene 
oxide and its application for real-time monitoring of protease activity w. ChemComm, 
pages 2399-2401, 2011. 


Rochefort, A. and Wuest, J.D. Interaction of Substituted Aromatic Compounds with 
Graphene. Langmuir, 1:210-215, 2009. 


Demircioglu, Z. and Kastas, A. and Buyukgungor, O. Theoretical analysis (NBO, 
NPA, Mulliken Population Method) and molecular orbital studies (hardness, che- 
mical potential, electrophilicity and Fukui function analysis) of (E)-2-((4-hydroxy- 
2-methylphenylimino)methyl)-3-methoxypheno. Journal of Molecular Structure, 
1091:183-195, 2015. 


Selvaraj, S. and Rajkumar, P. and Kesavan, M. and Gunasekaran, S. and Kumaresan, $. 
Experimental and theoretical investigations on spectroscopic properties of tropicamide. 
Journal of Molecular Structure, 1173:52—62, 2018. 


Mallineni, S. and Shannahan, J. and Raghavendra, A. and Rao, A. and Brown, J. and 
Podila, R. Biomolecular Interactions and Biological Responses of Emerging Two- 
Dimensional Materials and Aromatic Amino Acid Complexes. ACS Applied Materials 
£ Interfaces, 8(26):16604-16611, 2016. 


Prasongkit, J. and Grigoriev, A. and Pathak, B. and Ahuja, R. and Sheicher, R.H. 
Theoretical Study of Electronic Transport through DNA Nucleotides in a Double- 
Functionalized Graphene Nanogap. J. Phys. Chem. C., 117:15421-15428, 2013. 


107 


BIBLIOGRAFÍA 


[90] 


[91] 


[92] 


[93] 


[94] 


[95] 


[96] 


[97] 


[98] 


Souza, L. and Amorim, R.G. and Scopel, W. and Ralph.H Scheicher. Electrical detec- 
tion of nucleotides via nanopore in hybrid graphene/hh-BN sheet. Nanoscale, 9:2207— 


2212, 2017. 


Reichert, J. and Ochs, R. and Beckmann, D. and Weber, H.B. and Mayor, M. and Loh- 
neysenm H.V. Driving current through single organic molecules. Phys. Rev. Lett., 
88(17):176804—1, 2002. 


Yagmurcukardes, N. and Aydin, H. and Can, M. and Yanilmaz, A. and Mermer, O. 
and Okur, S. and Selameta, Y. Effect of aromatic sams molecules on graphene/sili- 
con schottky diode performance. ECS Journal of Solid State Science and Technology, 
5(7):69—73, 2016. 


Rodríguez, S.J. and Albanesi, E.A. Electronic transport in a graphene single layer: 
application in amino acid sensing. Phys.Chem.Chem.Phys., 21:597-606, 2019. 


Mao, $. and Lu, G. and Yu, K. and Bo, Z. Specific protein detection using thermally 
reduced graphene oxide sheet decorated with gold nanoparticle antibody conjugate. 
Adv. Mater, 22(1):3521-35326, 2010. 


Huang, Y. and Dong, X. and Liu, Y. and Lic, L.J. and Chen, P. Graphene-based 
biosensors for detection of bacteria and their metabolic activities. J. Mater. Chem, 
21(1):12358-12362, 2011. 


Dong, X. and Fu, D. and Xu, Y. and Wei, J. and Shi, Y. and Chen, P. and Li, J. Label- 
free electronic detection of dna using simple doble-walled carbon nanotube resitors. J. 
Phys. Chem, 112(1):9891-989, 2008. 


Hwang, M.T. and Landon, P.B. and Lee, J. and Choi, D. and Mo, A.H. and Glinsky, G. 
and La, R. Highly specific snp detection using 2d graphene electronics and dna strand 
displacement. PNA, 113(1):201603753, 2016. 


Shaibley, J.R and Yu, H. and Clark, G. and Rivera, P. and Ross, J. and Seyler, K. and 
Yao, W. and Xu, X. Valleytronics in 2d materials. Nature Reviews Materials, 1:16055, 
2016. 


108 


BIBLIOGRAFÍA 


[99] Chen, F. and Choun, M. and Chen, Y.R. and Wu, Y.S. Theory of valley-dependent 
transport in graphene-based lateral quantum structures. Condensed Matter and Mate- 
rials Physics, 94:112-119, 2016. 


109 


